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Abstract. We prove an asymptotic formula for the probability that, if one 
chooses a domino tiling of a large Aztec diamond at random according to the 
uniform distribution on such tilings, the tiling will contain a domino covering 
a given pair of adjacent lattice squares. This formula quantifies the effect of 
the diamond's boundary conditions on the behavior of typical tilings; in ad- 
dition, it yields a new proof of the arctic circle theorem of Jockusch, Propp, 
and Shor. Our approach is to use the saddle point method to estimate cer- 
tain weighted sums of squares of Krawtchouk polynomials (whose relevance 
to domino tilings is demonstrated elsewhere), and to combine these estimates 
with some exponential sum bounds to deduce our final result. This approach 
generalizes straightforwardly to the case in which the probability distribution 
on the set of tilings incorporates bias favoring horizontal over vertical tiles or 
vice versa. We also prove a fairly general large deviation estimate for domino 
tilings of simply-connected planar regions that implies that some of our results 
on Aztec diamonds apply to many other similar regions as well. 



1. Introduction 

1.1. Statement of the main theorem. Random domino tilings of finite regions 
often exhibit surprising statistical heterogeneity. Such heterogeneity would be ex- 
pected in the vicinity of the boundary, but in fact the presence of a boundary can 
make its influence felt well into the interior of the region. The research that led 
to this article is part of an ongoing effort to understand this phenomenon. The 
results proved here are the first to give a precise description of how local statistics 
for domino tilings can vary continuously throughout a region in response to the 
imposition of specific boundary conditions. 

Those who study random tilings of finite regions (in the plane) by dominos have 
tended to focus on regions that are rectangles of even area. In particular, Burton 
and Pemantle |BP[ | have done an intensive analysis of the small-scale structure 
of such tilings. Their work shows that once one gets away from the boundary 
of the rectangle, random tilings tend to exhibit statistical isotropy. Among all 
random processes that take their values in the set of domino tilings of the plane, 
the Burton-Pemantle process has maximal entropy, and it is unique in this regard; 
for this reason alone, it is worth further study. 
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Figure 1. A random domino tiling of an Aztec diamond of order 64. 



However, if one looks at random domino tilings of tileable finite regions in gen- 
eral, one finds that local behavior far from the boundary need not be governed by 
maximal entropy statistics, but can look very different. Moreover, the local behav- 
ior seen in one part of the region is in general different from local behaviors seen 
elsewhere. 

One especially tractable proving ground for the study of this statistical hetero- 
geneity has been the family of finite regions known as Aztec diamonds, introduced 
and studied in | EKLP | . Figure |l| shows an Aztec diamond of order 64 tiled randomly 
by dominos. In general, the Aztec diamond of order n can be defined as the union 
of those lattice squ ares whose interiors lie inside the region \{x, y) x + y < n + 1} . 

It was shown in | JPS| (and will be proved in subsection 6.4 by different methods) 
that, asymptotically, the circle inscribed in the Aztec diamond of order n serves 
as a boundary between domains of qualitatively different behavior. We call this 
circle the arctic circle, because, as one can see from Figure 0, the dominos outside 
the arctic circle are frozen into a brickwork pattern. To state the theorem more 
precisely, we impose a checkerboard coloring on the Aztec diamond of order n, so 
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that the leftmost square in each row in the top half of the diamond is white. We say 
a horizontal domino is north-going or south-going according to whether its leftmost 
square is white or black, and we say a vertical domino is west-going or east-going 
according to whether its upper square is white or black. (The motivation for this 
terminology comes from the "domino shuffling" algorithm introduced in | EKLP| 
and used in both | JPS| and |GIP|; this algorithm permits one to generate random 
domino tilings of Aztec diamonds in such a way that every possible tiling has the 
same probability of arising as every other, and indeed it was this algorithm that we 
used to generate the tiling shown in Figure |l|.) 

Say that two dominos are adjacent if they share an edge (i.e., their boundaries 
overlap on a segment of length 1 or more), and say that a domino is adjacent to 
the boundary of the Aztec diamond if it shares an edge with the boundary. We 
define the north polar region as the union of those north-going dominos that are 
each connected to the boundary by a sequence of adjacent north-going dominos. 
The south, west, and east polar regions are defined similarly, and the temperate 
zone is the union of those domin os th at belong to none of the four polar regions. 

The arctic circle theorem of JPS | states that for every e > 0, if one takes n 
sufficiently large, then for all but an e fraction of the domino tilings of the diamond 
of order n, the border of the temperate zone stays within distance en of the circle of 
radius n/\/2 with center (0,0). In particular, this implies that if one increases the 
radius of the circle by en, then with probability greater than 1 — e, in each of the 
four regions in the Aztec diamond that lie outside the enlarged disk, all dominos are 
aligned with their neighbors in brickwork patterns. The theorem also implies that 
if one decreases the radius of the disk by en, then within the shrunken disk dominos 
with different orientations are in some sense interspersed among one another (with 
probability greater than 1 — e); however, the theorem by itself gives no information 
on their distribution. 

In Theorem ^ of this article wc will give a quantitative analysis of the behavior of 
random tilings in the inner, disorderly zone. In particular, we will give an asymp- 
totic formula for the proportion of domino tilings of the Aztec diamond of order n 
that contain a domino at a specified location, i.e., the placement probability ioi that 
location. This formula depends only on the orientation of the domino, its parity 
relative to the natural checkerboard coloring of the Aztec diamond, and the relative 
position of the domino within the Aztec diamond (in normalized coordinates). One 
consequence of our formula is that random domino tilings exhibit "total statistical 
heterogeneity" within the central zone. That is to say, any two patches within the 
temperate zone that are macroscopically separated (i.e., separated by a distance 
on the order of n) will exhibit distinct statistics. (For a precise statement, see 
subsection 6.5.) 

Our work builds on the generating functions derived in GIP|. One of them is a 
rational function in three variables whose coefficients are the placement probabil- 
ities for which an asymptotic formula is sought. The authors of the earlier article 
carried out a relatively straightforward complex integration to calculate coefficients 
corresponding to dominos in the 2x2 block in the middle of the Aztec diamond; 
the resulting exact formula implies that in a diamond of order n, these placement 
probabilities are j + 0(i). In the present article we will apply the saddle point 
method to estimate contour integrals associated with more general coefficients of a 
related generating function (also derived in | ]GIP | ) . 
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We can now prepare to state our main result. We call the union of two adjacent 
squares in the Aztec diamond a domino space, to avoid confusion between actual 
dominos occurring in a particular tiling and the locations in which dominos can 
occur. Domino spaces are classified as north-going, south-going, west-going, or 
east-going in the obvious way, so that for instance a domino is north-going if and 
only if it occupies a north-going domino space. Because of symmetry, we lose no 
generality by focusing on the placement probabilities associated with north-going 
domino spaces. The midpoint of the bottom edge of each north-going domino space 
is some point {£, m) with \£\ + \m\ < n — 1. We call this the location of the north- 
going domino space. Normalizing by dividing by n, we obtain some point (a;, y) 
with I a; I -I- \y\ < 1. We call this the normalized location of the north-going domino 
space. 

Theorem 1. Let U be an open set containing the points (±i, i). // {x,y) is the 
normalized location of a north-going domino space in the Aztec diamond of order 
n, and (x, y) ^ U , then, as n oo, the placement probability at {x, y) is within 
o(l) ofV{x,y), where 

if x'^ + y'^ >\ and y < \, 
if x'^ + y"^ > i and y > ^, and 

ifx^ + < i. 

The o(l) error bound is uniform in {x,y) (for {x,y) ^U). 

Similarly, the south-going, east-going, and west-going placement probabilities 
near {x,y) are approximated by 'P{—x,—y), V{—y,x), and 'P{y,—x), respectively. 
This follows from Theorem |l| by rotational symmetry. 

The organization of the rest of this article is as follows. 

In the remainder of Section ^ we discuss some qualitative features of the main 
theorem and give some preliminaries for the proof. In Section we use the saddle 
point method to derive asymptotic estimates for certain numbers known as creation 
rates, which give placement probabilities when summed appropriately. In Section ||, 
we use this result to estimate, modulo an error term, the north-going placement 
probabilities. In Section ^, we use techniques from the theory of exponential sums 
to justify our bound for the error term. This completes the proof of the theorem 
away from the boundary of the diamond; Section |^ provides the final arguments 
that handle locations near the boundary. 

Section ^ discusses some consequences of the theorem. In particular, by taking 
a detour through the theory of domino tilings in general, we show that some con- 
sequences of the arctangent formula apply not only to the particular shape we call 
the Aztec diamond but also to slightly deformed versions of this shape (Proposi- 
tion ^0|) . We also give a new proof of the arctic circle theorem and a large deviation 
estimate for certain properties of random tilings of simply-connected finite regions 
(Theorem and Proposition ^) . Section ^ briefly sketches how the method of 
proof of Theorem |l| can be adapted to handle the more general case of random 
domino tilings when there is a bias in favor of one domino orientation over the 
other (horizontal versus vertical) . We conclude in Section ^ with speculations and 
open questions. 

For a treatment of the probabilistic preliminaries needed for Section @, see H . 



V{x,y) = 




tan 



2y-l 
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1.2. Features of the result. As a first comraent on the qualitative features of this 
formula, we point out the continuity of the formula for V{x, y) (except at (±i, i)). 
Indeed, if we had been so naive as to ask for an asymptotic formula for the placement 
probabilities for all horizontal domino spaces in an asymptotically small patch of 
the Aztec diamond (south-going as well as north-going), we would not get a single 
value at all but rather a pair of values, namely V{x,y) and V{—x,—y), which 
are not in general equal. That is, the local statistics are not even approximately 
invariant under translations that exchange the two color-classes. It is therefore 
all the more pleasant that the local statistics are asymptotically invariant under 
translations that preserve the two color-classes (at least, they are invariant if, in 
discussing local statistics, we confine ourselves to placement probabilities, and do 
not inquire about correlations between placements). 

Another important feature of the formula is the singular behavior that occurs 
near the normalized locations (±i, i), which we can explain as follows. In | EKLP| 
it is shown that the Aztec diamond of order n has exactly 2"("-+i)/2 domino tilings, 
and a formula derived in that article (formula (7) of Section 4) can be used to show 
that for < fc < n, exactly (^^^2"'"^^^/^ of the tilings have horizontal dominos 
covering the leftmost squares in the first k rows from the top and have vertical 
dominos covering the leftmost squares in the next n — k rows. Thus, the placement 
probability associated with the leftmost north-going domino space in the A;th row 
is exactly the sum 




This sum is very close to 1 for fc — ^ ^ —y/n and very close to for fc — ^ ^ y/n; 
macroscopically speaking, the placement probability jumps from 1 to discontin- 
uously. It might be possible to analyze the limiting behavior of the placement 
probabilities in the vicinity of the singularities under suitable scaling, but we do 
not explore this possibility here. 

An easily-understood symmetry property oi V{-,-) is the fact that 

(1.1) rix,y) = Vi-x,y). 

This is a consequence of the fact that reflecting a domino tiling through the line 
a; = carries north-going domino spaces to north-going domino spaces. A further 
identity satisfied by 7^(-, •) is the relation 

(1.2) V{x, y) + V{~y, x) + V{^x, -y) + V[y, -x) = 1. 

To see why this is true, one need only observe that the four domino spaces that 
contain a particular lattice square (fewer, if the square is on the boundary) must 
have placement probabilities that sum to 1. 

A subtler consequence of Theorem ^ is the fact that the level sets of 7^(x, y) (for 
probabilities strictly between and 1) are arcs of ellipses. More specifically, for < 
p < 1 the level set {[x, y) : V{x, y) = p} and the level set {(x, y) : 7'(x, ?/) = 1 — p}, 
together with the singular points (±5, 5), form an ellipse tangent to the boundary 
of the diamond at the two singular points. As p ^ (or p 1), the ellipse becomes 
the inscribed circle, which is the zero-set of the function 2x'^ + 2y^ — 1; in the case 
p = \, the ellipse degenerates into the line segment joining the two singular points, 
which is the part of the zero-set of the function (2y — 1)^ lying inside the Aztec 
diamond; and in general, the ellipse will be the zero-set of some convex combination 
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Figure 2. Level curves of north-going placement probabilities. 



of + 2%)"^ — 1 and {2y — 1)^. The point (0, 0) hes on the level set p = j, which is 

3 

4' 



an arc of an ellipse; the complementary arc of the ellipse is the level set p — |, and 



the point on this arc opposite (0,0) is the point (0, |). The situation is depicted 
schematically in Figure |^. 

1.3. Preparation for the proof. Recall that, under the original (unnormalized) 
coordinate system, each north-going domino space in an Aztec diamond of order 
n is assigned some location {£,m) with \£\ + \m\ < n — 1. It is easy to check that 
i + m must have the same parity as n — 1. Define P(^, m;n) as the probability 
that a random domino tiling of the Aztec diamond of order n will have a domino 
occupying the north-going domino space at location (£, m); for \£\ + \m\ > rt — 1, or 
£ + m ^ n — 1 (mod 2), define P{£,m;ri) = 0. For instance, we have P(0,0; 1) — ^, 
P(0, 1; 2) = |, and P(0, -1; 2) = P(l, 0; 2) = P(-l, 0; 2) = i. 
Define 

(1.3) Ct{£, m; n) = 2{P{£, m; n) - F{£, m - 1; n - 1)). 

This quantity is called the net creation rate at location (£,m), but the reason for 
this name and the interpretation of the quantity in terms of domino shuffling are 



not needed for our purposes. (For the motivation, see |GIP|.) 

Define c(a, 5;n) to be the coefficient of z°- in (1 -I- z)"~^(l — z)''. (Note that 
c(a, b; n) is the Krawtchouk polynomial Pa evaluated at 6. For information about 



LOCAL STATISTICS FOR RANDOM DOMINO TILINGS 



7 



Krawtchouk polynomials, see [MS, p. 130].) Our proof of Theorem |T] will be based 



on the following result from |GIP| 



Proposition 2. Let n> 0. Suppose £ and m are integers with £ + m — n (mod 2) 

and \£\ + \m\ < n. If we let a = {£ + m + n) /2 and b = {£ — m + n) /2, then 

Ct{£, ■m;n + l) = c{a,b; n)c{b, a; n) /2". 
For other integers £ and m, we have Ct{£, m; n + 1) = 0. 

This proposition implies that the creation rates are non-negative, if we use the 
identity c(6, a;n)6!(n — 6)! — c{a,b;n)a\(n — a)!. (This identity is a standard 
fact about Kr awtc houk polynomials, and follows immediately from Theorem 17 
on page 152 of |MS|.) When combined with Proposition H, the identity implies that 



Cr(£, TO;ri + 1) is a positive factor times the square of a Krawtchouk polynomial 



and hence that Cr(£, m;n + 1) > 0. Note that this inequality, together with (1.3), 
yields 

(1.4) P{£,m]n) <P{£,m + h;n + h) 

for ft. > by induction on h. 

We will also need the following result on exponential sums. 

Theorem 3 (Kusmin-Landau). Let \ \ ■ \ \ denote the distance to the nearest inte- 
ger, I be an interval, and f be a real-valued function on I . If f is continuously 
dijjerentiable, f is monotonic, and ||/'|| > A > on /, then 

exp{27Tif{n)) = 0{X-'). 

nG/nZ 

The constant implicit in the 0{X~^) term does not depend on I. 



A proof can be found in [GK, p. 7] 



1.4. Out line of the Proof of Theorem 1. We begin the proof of Theorem |^ 



by using (O) to write the placement probabilities as sums of creation rates, which 



gives the formula 

(1.5) P(^,m;n) = i^Cr(£,TO- fc). 

fe>0 

(Note that the remark after Proposition |^ shows that, as claimed in the abstract, 
this is a weighted sum of squares of Krawtchouk polynomials.) We will estimate 
the creation rates, and then use our estimate to prove the asymptotic formula for 
the placement probabilities. 

Because of Proposition ^, to estimate the creation rates it suffices to approxi- 
mate the coefficients of the polynomials (1 -I- z)"~''(l — z^. To do this, we write 
the coefficients as contour integrals in the usual way, and then apply the sad- 
dle point method to these integrals. Sufficiently far outside the arctic circle, this 
method shows that the creation rates are exponentially small in n (Proposition ||) ; 
sufficiently far inside, it approximates them by a well-behaved function times an 
oscillating factor (Proposition ^) . With additional work, it might be possible to 
obtain a uniform estimate over the entire Aztec diamond, but we can make do with 
just these estimates. 



We would then like to substitute our creation rate estimates into (1.5) and con- 



vert the sum to an integral to determine its asymptotics. If we are willing to be 
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unrigorous, we can wishfully replace the oscillatory cosine-squared factor in Propo- 
sition ^ by its mean value i , obtaining (for locations inside the inscribed circle) 



P{i,m;n) 



k=0 




fc)2 - 2^2 „ 2(m- fc)2 

k + n — 2m 
- fc)2 - 2^2 _ 2{m~ky 
n — 2m 



dk 

k=0 



\/ — 2P — rn? 

2y-l 
v/1 - 2a;2 - 2y2 ' 

where k — tmax is the larger of the two roots of the equation (n — fc)2 — 2£2 _ 2(m — 
fc)^ = 0, a; = ^/n, and y = m/n] we truncate the sum and integral at t,nax on 
the supposition (to be discussed in the next paragraph) that essentially no creation 
occurs outside the arctic circle. To make this argument rigorous, we need to deal 
honestly with the oscillating factor in the creation rate estimate inside the arctic 
circle, and we need to circumvent the non-uniformity of our estimates. 

The non-uniformity of the estimates can be dealt with simply by summing over 



a smaller interval than in (1.5). The exponentially small bounds on the creation 



rates outside the arctic circle show that the terms in (1.5) that come from locations 
outside the arctic circle contribute very little to the sum. Motivated by this, we 
look at the sum of all the terms that come from locations that are far enough inside 
the arctic circle that our creation rate estimates from Proposition ^ apply. Because 
the creation rates are all non-negative, this new sum underestimates the placement 
probability. Dealing appropriately with the oscillating factor (as described below) 
gives an estimate for the new sum; as we see in the computation above, this estimate 
turns out to be the arctangent formula from Theorem |l|. A short argument shows 
that the placement probabilities can be no larger asymptotically (Proposition |l2|) , 
and because our estimate is an underestimate we know they can be no smaller. 
This completes the proof. (Actually, this method works only away from the edges 
of the diamond, so it is not until Section |^ that the proof is completed.) 

All that remains is to describe how to deal with the oscillating factor in the 
summand. We must show that replacing the oscillating factor by its average value 
has an asymptotically negligible effect on the sum. Equivalently, we must show 
that the difference between the original sum and the smoothed sum is small. This 
difference is an exponential sum, and we can estimate it using the Kusmin-Landau 
Theorem once some preparatory results (Lemmas ^, 0, and ^4|) are in place. 

2. Creation Rate Estimates 

Our proof of the asymptotic formula for placement probabilities begins with an 
estimate of creation rates, which is proved using the saddle point method. Because 
Proposition |^ is most conveniently stated for an Aztec diamond of order n + 1 , we 
will estimate the creation rates in an Aztec diamond of order n + \. From this 
point on, we assume that i + m = n (mod 2), because otherwise Cr(£, m; n + 1) is 
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necessarily 0. As pointed out in subsection 1.4, the creation rates behave differently 
inside and outside the inscribed circle. If we estimate the creation rates inside it, 
we get the following result: 

Proposition 4. Fix e > 0. If i'^ + m'^ < {1 — e)n'^ /2 and ( + m = n (mod 2), then 

for some function ^{£,m;n), which we determine explicitly below. 

The subscript in Oe{n~'^) indicates that the implicit constant depends on e. In 
this paper, if any subscripts appear on a big O term, then the implicit constant 
depends only on the indicated variables, but the absence of subscripts should not 
be taken to imply that the implicit constant is absolute. 

Proof. Let 

(l + z)"-''(l-z)'' 



(2.1) f{z) 

where a and b are defined as in the statement of Proposition ||. To approximate 
the creation rate, we need to approximate c(a, 6;n), which is the constant term 
of f{z). The constant term is given by the usual contour integral, which we will 
approximate using the saddle point method. 

Write a = (1 +u)n/2 and 6 = (1 + v)n/2, so that —l<u,v<l. Note that the 
u, V coordinates are related to the coordinates in the statement of the proposition 
by u = + m) /n = X + y and v = {£ — m) /n = x — y. Wc will keep u and v fixed 
as we send n to infinity. 

The critical points of f{z) are 



and 

Z2 



1 - u 



+ t;^ — 1 



1 - u 

Because ^2 + 771^ < 71^/2, we have u^ + w^ < 1. It follows that zi and Z2 are complex 
conjugates on the circle |zp = (1 + u)/{\ — u). We now apply the saddle point 
method. To find the constant term of f{z), we integrate f{z)/{2TTiz) about the 
circle of radius a/(1 + u)/(l — u) centered at the origin. One can check that, on this 
circle, \f{z)\ is greatest at the critical points zi and Z2. (To check it, parametrize 
the circle by the angle 9 formed with the real axis. One has d\og\f{z)\^ /d9 = 
iff z is one of the two critical points or z lies on the real axis. At the critical 
points, 92 log |/(z)p/902 = n(u2 + - 1)/(1 — w^) < 0, so \ f{z)\ has maxima at 
these points. It must have minima on the real axis, since between any two local 
maxima there must be a local minimum.) As n goes to infinity, the integral is given 
asymptotically by the integrals over the parts of the path near the critical points, 
which can be estimated straightforwardly. This is the saddle point method. We 
will omit the details of the argument leading to the approximation, because they 



are standard, and can be found, for example, in [dB, pp. 87- 
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The saddle point method teUs us that the constant term of f{z) is the sum 
Zi{l + 0(n"i)) + Z2{1 + O(n-i)), where 



(2.2) 

and 

(2.3) 



Z2 — Z\ 



/(fi) / Stt 
2^zi Y (log/)"(zi) 



27r 



27rz2 V (log/)"(;22) 



(For the proof of Proposition ^ we wiU not need to determine the signs of the square 
roots in (2.2) and (2.3), but they must be chosen so that Z^ and Z-i are complex 
conjugates.) 

Simplifying z^(log/)"(z) yields 



n[l + U'~4:Z^-2 uz^ - 4 vz^ + uz^ - z^) 



2(z2_i)2 

From this, one can check that at either critical point of f{z), z'^{log f)" (z) has 
absolute value 



(2.4) 



\zf {log fnz.)\ = 



n / (1 — — u2)(l _ u2~) 



(l~f2) 



Let ^{u,v;n) be the phase of Zi, so that Zi — |Zi| exp(i ^(u, w; n)). Then 
Zi + Z2 = 2|Zi|cos*(w,u;n), 
and c(a, 6; n) is approximated by 

\Zi\ 



(2.5) 
Of course, 

(2.6) 



c{a,b;n) = 2\Zi\ cos^ {u,v;n) + O 



\ZA = 



27r 



l/(^i)l / 

2n y \zf {log fY'{z,)\- 



Since |1 + zip = 2(1 - v)/{l ~ u), |1 - zip 2(1 + v)/{l - u), and |zip 
(1 + u)/{l — u), we see that 



(2.7) 



\f{zi)\=2-/' 



(l„-y)("-ft)/2(l+^,)fc/2 



(1 - m)(«-'')/2(1 +w)''/2' 

We have Cr(£, m; n + 1) = c(a, 6; n)c{h, a; n) /2", by Proposition ||. Interchanging 
a and h corresponds to interchanging u and v. Let i denote the result of inter- 
changing u and V (and also a and 6) in the expression x, so that, for example, 



u — 2v = V — 2u. When we substitute ( |2.7| ) and ( p.4| ) into (2.(:), we see that 

, on 

l^lll^ll- , 

7mv 1 — u — w 

Hence, by (U) 

4cos5'('u,t;;n)cos^'(w,u;n) _2. 

Cr(£,TO;n+l) = , 1- 0{n ). 

Trnv 1 — — 
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(To see that the error term is 0{n one uses the fact that it is 0(|Zi||Zi|/(n2")) 
and that |Zi||Z^| = 0(2"/n).) 

Now we check that cos ^'(w, u; n) = ± cos v; n). The identity 

c{b, a; n)b\{n — 6)! = c(a, b; n)a\{n — a)! 

suggests that this should be so, but does not seem to prove it. If we set a = 
zf{\ogf)"{zi), we find (after some simpUfication) that 

a 1 — 



a 1 — 

Thus, a and a have the same phase. If we combine the formulas 

1 + Zl _ 1-v 
1 + ^1 1 — M 

and 

1 + Zl 1 — V Zl 

with 

f{zi) = {i+zir' 



1 + Zl J \Zi 



we find that f{zi) equals /(zi) times a positive factor, so their phases are equal. 
Because 

Zi^±lM,[^ and Zi^±l^^ 



2tt y a 2tt 

we see that Zi and Zi have the same phase, to within a sign, so cos^(u,u;n) — 
± cos ^'(w, v; n). 

Finally, we change to the coordinates of our generating function by the substitu- 
tions u = (£ + 'm)/n and v — {£ — m)/n. We set $(^, m; n) — ^'(u, w; n). Then when 
I -\- m = n (mod 2), the creation rate at the (£, m) location in an Aztec diamond 
of order n + 1 is 

r (0 ±4cos2$(Vn^^„, _2^ 

Crli:, m; n + 1) = — + (J(n ). 

Because creation rates must be non-negative, the ± sign in this formula can always 
be taken to be 

The constant implicit in the big O depends continuously on u and v. Thus, for 
fixed e > the constant can be chosen uniformly for all u and v with u^ + f^ < 1 — e. 
We have therefore proved the result claimed in the statement of the proposition. □ 

Given m, and n, define 

(2.8) 5^ = {fc e Z : /c > and ^2 _^ (m - kf < (1 - e)(n - hfjl}. 

Also, define fcmax to be the greatest element of 5^, and fcmin to be the least. (In 
Section H, we will sum the creation rates in ( p_.5| ) as fc varies over Se- We will do so 
to make it possible to apply Proposition to the terms in the sum, as described in 



subsection 1.4.) 
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Lemma 5. Suppose that \£\ + |r7i| < (1 — 6)n for some fixed 6 > 0. Then fcmax < 
(1 — (2 — V^)S + 0{e))n. Hence, if e is small enough compared to 5, then [n — 
^max) "'^ = 0{n^^), and for k we have 

0(e-i/2n-i) = 0,,5(n-i). 



Proof. We have fcmax = (2y — 1 + a/2((1 — yY — x"^) + 0{e))n, where x = £/n and 
y — m/n. For fixed y, this function is clearly maximized at a; = 0. When a; = 0, 
it becomes a linear function maximized at y = 1 — 5 (for |2;| + |y| < 1 — S). This 
yields fcmax < (1 - (2 - V2)S + 0{e))n. Therefore, (n - fcmax)"^ = 0(n"i) if e is 
small enough compared to 6. 
For k G Se, we have 



It follows that 



Therefore, 



£2 + (m - fc)2 <{n- fc)2(l - e)/2. 
{n - fc)2 - 2e^ - 2(to - fc)2 >in- kfe. 



^{n-kf -2P -2{m-kY 



If we are not worrying about dependence on e, this is 0(n ^). This completes the 
proof. □ 

In order to apply the Kusmin-Landau Theorem to deal with the exponential 
sums that will appear later in the proof (as discussed in subsection 1.4), we will 
need to specify $, since the phase of Zi is not uniquely determined. Also, it will be 
convenient to extend it to a function of real, and even complex, variables (rather 
than just integers). 

Given a point {x, y) ^ (0, 0) in the plane, define 0{x^ y) to be the angle in (— tt, tt] 
formed by the right half of the horizontal axis and the ray from the origin through 
{x,y)- 

Lemma 6. We can choose $(£, m; n) in Proposition so that if one sets I = xn, 
m = yn, and k = nn in d^{£, m — k; n — k)/dk, then d^{i, m — k; n — k) / dk equals 



e 



(^-x + y - K, v/(l - k)2 - 2 {x^ + {y - k)2)^ - 

e(l-K-2x, - k)2 - 2(.t2 + [y-nf)^ + 
a;2 — K — 3 j/K + 2 k2 + + 2/2 
n^(l - k)2 2 (a;2 + (y - k)2) (y + l- 2K-a;)(2/ + l- 2K + a;)' 

As n tends to infinity, the last term is 0{l/n) for fc € 5^ with e > fixed, as long 
as \i\ + \m\ < (1 — S)n for some fixed S > 0, and e is small enough compared to S. 



Proof. From ( |2.2[ ), we see that we can choose $(£, m; n) to be the imaginary part 
(2.9) Im ( log f{zi)- log zi - ilog((log/)"(zi)) 
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If we substitute m — k for m and n — k for n and differentiate, then the first term 
of (2.9) contributes the ^-terras in the formula we are proving. To see this, recaU 
that (up to an irrelevant multiple of 2ni) 

log/(zi) = (n - 6)log(l + zi) + 61og(l - zi) - alogzi. 

After we express this in terms of n, £, and m and substitute m — k for m and n — k 
for n, the right hand side becomes 

{n/2 -k-e/2 + m/2) log(l + zi) + {n/2 + 1/2 - m/2) log(l - ii) 
(2.10) -(n/2-fc + ^/2 + m/2)logzi, 

where ii is the function of k that results from making the substitutions in zi. 
Denote by L the function ( 2.1C ) . When we differentiate L with respect to k (holding 
n, £, and m fixed), we get 

, . , . , dL dzi 

- = log.i-log(l + zi) + ^^. 

Because Zi is a critical point of /, dL/dzi = 0, so dL/dk — logzi — log(l + zi). 
Now expressing the imaginary parts of the logarithms in terms of 9 gives the desired 
terms from the formula we are proving. (To simplify the terms to the form found 
in the statement of the lemma, one has to use the fact that for a > 0, 9{ax, ay) = 
0{x,y).) 

When we substitute and differentiate, the remaining terms in (2.9) clearly give 
algebraic results. We omit the details of the calculations, since they are tedious 
and straightforward. 

The claim that the last term is 0{l/n) for fc e S'e is a consequence of Lemma |[ 
The only thing to check is that although the denominator vanishes at k = (y + 1 ± 
a;)/2, these two points are nev er in (or near enough to cause problems). To see 



that, note that the definition (2^) oi is equivalent to the set of fc > for which 



(2.11) (1 - k)2 ^ 2{x^ + iy- k)2) > £(1 - k)2. 

Note that k = 1 is impossible (since then we must have x — and y — 1, so 
\x\ + \y\ > 1 — S). However, substituting k = {y + 1 zL x)/2 in the left hand side of 
( ^.11| ) gives -(3a; ± (1 - 2/))^/4. Thus, the factors y+l-2K-x and y+l-2K + x 
in the denominator of the last term in our main formula cannot become arbitrarily 
small, and the last term is indeed 0(l/n). □ 

In Section ^, we will need the following result (to make i t po ssible to apply 
exponential sum techniques in the way described in subsection 1.4). 

Lemma 7. The function (P^{£,m — k;n — k) / dk^ is algebraic. For any fixed n, I, 
and m satisfying \£\ + |m| < n and e > 0, there exists a neighborhood U in C of the 
smallest real interval containing such that as a function of k, d^^{£, m — k;n — 
k)/dk'^ is holomorphic on U . 

Proof. We will use the formula for d^{£, m — k; n—k)/dk from Lemma |^. Let U he a 
small, simply-connected neighborhood in C of the smallest real interval containing 
Se, such that the points k = n{y + 1 ± a;)/2 are not in U. (We checked at the 
end of the proof of Lemma || that these points are not in S^.) It follows from the 
definition of that 

(n - kf - 2{f + (m - kf) > e{n - kf > 
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on S'j. If n ~ k, then £ — and m ^ k — n, so \£\ + \m\ = n (contradicting 
1^1 + |m| < n). Thus, (n — k)^ — 2(£^ + (m — k)^) > e on S^, and hence there is a 
holomorphic square root of {n — k)'^ — 2(£^ + (m — fc)^) on J7, if U was chosen to be 
sufficiently smah. It follows that the third term (the algebraic term) of the formula 
for d^{£, m — k; n — k) /dk in Lemma ^ is holomorphic on U. The derivative of that 
term is thus algebraic and holomorphic on U , so to complete the proof we just need 
to check this for the other two terms. 

The first two terms can be expressed in terms of the arctangent. If we do so, we 
find that the derivative with respect to k of the sum of those two terms is 

ny/{l - k)2 - 2(a;2 + (y - K)2)(y + 1 - 2k - x){y + 1 - 2k + x)' 

This is also algebraic and holomorphic on U. Thus, d^^{£,m — k;n — k)/dk'^ is 
holomorphic on U and algebraic, as desired. □ 

We now know everything we need to know about how the creation rates behave 
inside the arctic circle. Outside the arctic circle, we can get an exponentially small 
upper bound for the creation rates. This will be used for bounding the placement 
probabilities outside the arctic circle (Proposition ^3|) . 

Proposition 8. For each e > 0, there exists a positive constant r < 1 such that 
whenever £^ + rr? > (1 + e)n^ /2, 

Cr(^,m;n + 1) = Cl(r"). 

Proof. We assume that £ + m = n (mod 2), since otherwise Ci{£, m; n + 1) = 0. As 
in the proof of Propositi on ^ , we will integrate f{z)/{2iTiz) around a circle about 



the origin, where, as in (2.1), 

, _ (1 + zr-\\ - zf 
nz)- 

However, since we are looking only for an upper bound and not for an asymptotic 
estimate, we will not need the full saddle point method. We will only sketch the 
proof, because the details are straightforward but somewhat tedious to check. 

We will use the same notation as in the proof of Proposition ^; for example, we 
write a = (1 + u)n/2 and 6 = (1 + v)n/2. Since + v"^ > 1 + e, the critical points 



Zl 

and 

Z2 



of f{z) are real. (Of course, the case u = 1 has to be handled separately, but this 
will not cause problems.) We will integrate f{z)/{2TTiz) around a circle of radius 
R, where R will be either \zi\ or \z2\- We choose R = \zi\ where |/(zi)| is the lesser 
of |/(.2i)| and 1/(2:2)1. To bound the integral, we will use the fact that the absolute 
value of the integral is at most as large as the greatest value of \f{z)\ on the circle. 

It is not hard to check by straightforward manipulation of inequalities that 
\f{zi)\ > \f{z2)\ if uv > 0, and |/(zi)| < |/(z2)| if uv < 0. (Since \u\,\v\ < 1 
and + v"^ > 1 + e, we cannot have uv = 0.) Thus, R — \z2\ if uv > 0, and 
R— \zi\ otherwise. 



LOCAL STATISTICS FOR RANDOM DOMINO TILINGS 



15 



Take i G {1,2} so that R — \zi\. On the circle of radius R about 0, \f(z)\ is 
greatest when z — Zi; in fact, the second derivative test shows that this is the only 
local maximum. Thus, the integral is bounded by \ f{zi)\, so \c{a,b;n)\ < \ f{zi)\. 

Because the sign of uv doesn't change when u and v are interchanged, 1/(2:^)1 is 
the lesser of [/(^i)] and |/(z2)|. Hence, \c{b,a;n)\ < It follows that 

C.(, + MM. 

A simple calculation gives |/(2i)| 1/(22)! = 1/(^2)11/(^1)1 = 2". The inequalities 
|/(zi)| < \f{z3-i)\ and |/(^i)| < |/(2;3_i)|, together with the fact that the only 
dependence on n in any of these expressions is in the exponent, imply that the 
creation rate at {u,v) is 0(r") for some r < 1. A little more care in the estimates 
shows that this bound can be chosen uniformly for u'^ + v'^ > 1 + e, as desired. □ 



3. Placement Probability Estimates 

Now that we know the creation rates, we can determine the placement prob- 
abilities. Fix 6 > 0. In this section, we will look at the placement probabilities 
P(^, m;n-|-l) at points {£,m) satisfying £ + m = n (mod 2) and 1^1 + |m| < {1 — S)n. 
(The congruence condition rules out the placement probabilities that we know are 
0, and the inequality lets us apply results such as Lemmas 5-7.) 



From (1.5), we see that 

P(^, to; n + 1) = i ^ Cr(^, m - k;n + l~ k). 

k>0 

It will turn out that the creation rates on or beyond the arctic circle make a van- 
ishing contribution to this sum as n —s- 00, so we can remove them from the sum 
without affecting its asymptotics. To remove these terms from the sum, fix e > 
(which we assume is small compared to 6, so that we can apply results such as 
Lemma ^), and look at the sum 

Pe = ^ ^ Ci{£, m - k;n + l~ k), 



where is defined by (2.8). (Note that sometimes = 0; we will see that this 
occurs only when the placement probability is nearly 0.) We will approximate P^, 
and prove that it approximates P{£, m;n + 1). First, we prove a few easy lemmas. 

Lemma 9. Consider the equation (1 — t)^ — 2x^ — 2(y — t^ ~ 0. For \x\ + |?;| < 1, 
this equation has two real roots t. The greater root is iff + = 1/2 and 
y < 1/2, and is less than iff x^ +y'^ > 1/2 and y < 1/2. The lesser root is greater 
than or equal to iff + > 1 /2 and y > 1/2. 

Proof. Since the discriminant of the polynomial is 8{l — x~y){l + x — y), we see that 
it has two real roots whenever \x\ + \y\ < 1. Clearly, is a root iff x^ + y'^ = 1/2, 
and since the sum of the roots is 4?/ — 2, it is the greater root iff also y < 1/2. One 
can check the other claims similarly. □ 



Lemma 10. Let i5 > 0, and suppose \x\ + \y\ < I — 5 . Let k{s) be any branch of 
the multivalued function of e defined by (1 — £)(1 — k)^ — 2x^ — 2{y — k)^ = 0. Then 
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for £ > (and sufficiently small relative to 5), we have K{e) — k(0) + Os{£)- (The 
constant implicit in the Os{£) does not depend on x, y, or e.) 

Proof. This simply amounts to showing that ^'{e) is bounded as a function of x, 
y, and e, for e sufficiently small. If one computes K(e) using the quadratic formula, 
and then differentiates it with respect to e, one finds that it equals a continuous 
function of x, y, and e (for e near 0) divided by 



^{\-e){y-lY-{l + e)x\ 

If £ is small enough compared to (5, then K'(e) will be continuous, and hence 
bounded, for all x, y, and e with |a:;| + \y\ < 1 — ^. Then K(e) ~ k(0) + 0(e), 
as desired. □ 



Proposition 11. Let 5 > Q and e > 0, such that e is sufficiently small compared 
to 6. If \£\ + \m\ < (1 - S)n and l + m = n (mod 2), then 

P, = V{i/n, m/n) + 0,,s{n-^) + Os{e^'^). 

Proof. From formula ( |1.5| ) and Proposition |^, we see that P^ is approximated by 

O(n-i). 



1 4cos2 $(£, m - fc; n - fc) 



2 ny/{n - k)^ - 2P ~ 2{m - ky 

(The error here is bounded by a constant multiple of X!j>ri-fc J~^i which is 
0((n — fcniax)^^), and hence 0{n~^) by Lemma ||.) 

Since 4cos2z — exp(2iz) + exp{—2iz) + 2, we see that is given to within 
0{n^^) by the sum of 

1 ^ 2 

2 ^ ny/{n- kf -2P -2{m-W 



(3.1) 
with 



:exp{2i^{£,m— k;n— k)) 



2 ^V(^-fc)'-2f'-2(m~fc)2 

and its complex conjugate. Proposition |l^ will show that the latter two sums are 
0{n~^) as n goes to infinity. 

Assuming Proposition 16, we can prove the desired limit by approximating the 
sum (^^) with an integral. The sum is equal to 

(3.2) If— 2 :dk + 0{n-^); 

2jk^,„ 7r^(n-fc)2-2<'2-2(m-fc)2 

to see that the error is 0(n~^), note that the summand (viewed as a function of a 
real variable k) is 0{n~^) by Lemma ^, and is monotonic on H {—oo,2m — n) 
and n (2m — n, oo). 

By Lemma ^, the polynomial {n — i)2 — 2^2 _ 2{m — t)'^ has real roots t. Let 
^min be the lesser root, and i,„ax the greater root. Then C [0, imax] if imax > 0, 
and = if tmax < 0. By Lemma we have tmax = iff ^2 _|_ ^^i^ = n'^/2 and 
m < n/2, and tmax < iff i'2 + m'^ > n2/2 and m < n/2. In both of these cases, we 
have V{l/n,m/n) — and P^ = 0{n~^). Thus, we need only deal with the case 

^max ^ 0. 
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Suppose iniax > and ^min < 0, i.e., P + < n?/2. It follows from Lemma [lO| 
that fcmin — 0{ne), and fcmax = ^max + 0{ne). We will approximate the integral in 
(il)by 

(3.3) - / =dk. 

2 Jo -K^in-kY - 2P - 2{m-kY- 

This approximation introduces further error. To see how large the error is, first 
rescale by a factor of n, so that the function under the square root sign becomes 
(1 — k)2 — — 2(?/ — k)^. Around a root r, this function can be expanded as 
— (k — r)2 ± 2-^/2^/ (y — 1)^ — x^{k. — r) (with the sign depending on which root r 
is). Because |a;| + < 1 — ^, the coefficient of k — r cannot become arbitrarily 
small. Thus, for small e, the error introduced by the approximation is bounded by 
a constant (depending on S) times 

de' 



and hence by 0(£^/2) 



One can evaluate the integral (3.3) explicitly, because 



,„ , dk 1 / fc + n — 2to 

(3-4) / ^========= = tan' 



^{n - fc)2 - 2^2 _ 2(m - kf \ ^ {n - kf - 2P - 2(m - k^ ) 

As k imax, the right hand side of (^^) approaches ^ (since the numerator of the 
fraction is positive as its denominator approaches 0). We see that ( |3.3| ) evaluates 
to 

1 1 J / 2™ — n 



— tan 



2 TT \^n^ - 2^2 _ 2m^, 

The case with tj^in > (i.e., i'^ + m? > and m > n/2) is completely 

analogous, except the integral is over the interval [imin, ^max], rather than [0,t,„ax]- 
This integral is 1, so we get that = 1 + 0{n^^) + 0{e^''^), which agrees with 
V{t/n,m/n) = 1. This proves the desired result. □ 

We still need to prove that approximates P(£, to; n + 1). We do that as follows: 

Theorem 12. Let 5 > Q and e > 0, such that e is sufficiently small compared to 
6. If \i\ + \m\ < (1 - S)n and £ + m = n (mod 2), then 

P(£, to; n + 1) = ne/n, m/n) + Osis^^^) + O^.^sin-^)- 

Proof. We need to show that P^ approximates P{i, m;n + 1). Since Proposition || 
implies that the creation rates are all non-negative, and P^ is the sum of a subset 
of the creation rates appearing in the sum giving P(^, TO;n + 1), the placement 
probability must be at least P^. 

Also, given any point in the Aztec diamond, the north-going placement prob- 
abilities at the four points obtained by rotating it by multiples of 90° about the 
origin sum to 1 . This is true because by rotational symmetry these placement prob- 
abilities are equal to the placement probabilities in each of the four directions at 
the original point, which must sum to 1. This is the content of ( |l.2| ), except here 
it is expressed in terms of the placement probabilities, rather than the asymptotic 
formula. 

One can check by direct computation that Vix, y) + 'P{y, —x) -\- V{—x, —y) + 
V{~y^ a;) = 1. If the difference between Ve and the placement probability were not 
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0(e^/^) + 0{n then the four placement probabilities would have to sum to more 
than 1, which is impossible. □ 

The statement of Theorem |2| implies that away from the edges of the diamond, 
the placement probabilities converge uniformly. (Given any e > 0, the theorem 
implies that if n is sufficiently large, then the placement probabilities are within 
a constant multiple of e^/^ of their limiting values. In fact, the slightly awkward 
theorem statement is equivalent to asserting uniform convergence; we state it that 
way because it seems to be the form in which it is most naturally proved, given our 
setup.) Thus, assuming Proposition [l^, we have very nearly proved Theorem |l|. In 
Section ^ we will complete the proof, using the following proposition: 

Proposition 13. For each e > 0, there exists a positive constant r < 1 such that 
whenever + > (1 + e)n'^ /2, 



P{£,m;n + 1) 



0(r") if m < n/2, and 

l + 0(r") ifm>n/2. 



Proof. First suppose that m < n/2. The desired result will follow from the equation 
(3.5) P{e, m; n + 1) = i ^ Ct{£, m - fc; n + 1 - fc), 

fc>0 

together with the estimate given by Proposition^. First, we show that Proposition|^ 
applies to the creation rates appearing in the sum. Consider 

as a function of k. Its first derivative at is 

+ m? — mn 



2- 



n 

which is greater than since P + m? > 12 > mn. The only root of the derivative 
is 

+ m^ — mn 
< 0. 



3 



Thus, the function (3.6) is increasing for < fc < n. (Note that in (3.5) we need 
only sum up to fc = (m + n)/2^ since beyond that point m — k < —{n — k) and 
hence Cr(^, m — fc; n + 1 — fc) =0. Thus, fc never reaches the pole in ( |3.6[ ) at n.) 
Therefore, £'^ + {m — fc)^ > (1 + £)("- ^ fc)^/2, and Proposition || appUes to bound 
the creation rates in (3.5). 

Thus, for some constant s between and 1, 

(m+n)/2 

P(£,m;n+1)< ^ 0(s"-'=). 

fc=0 

This geometric series is bounded by 0(s^"^'"-'/^) = 0(s"/'*). This proves the desired 
bound, with r = s^/^. 

For m > n/2, we use the trick of summing the placement probabilities at the 
four points obtained by rotating by multiples of 90° about the origin. As in the 
proof of Theorem the sum must be 1, and we know that three of the terms are 
0(r"). Therefore, the fourth must be 1 + 0(r"), as desired. □ 
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4. Exponential Sums 

In the proof of Proposition |ll|, we needed to show that is within 0{n^^) 
; to do so, we made use of an estimate whose proof was deferred 
. In this section, we will derive that estimate. We begin with the 
following lemma. 

Lemma 14. Let F(xi, . . . ,Xn+i) be an algebraic function ofn+1 variables (defined 
on a subset o/C"^^ to be specified shortly), and let S be a subset ofC". Suppose that 
for each (yi, . . . ,?;„) G S, there exists an open set U d C such that as a function 
of Xn+i, F{yi, . . . ,yn,Xn+i) is (defined and) holomorphic on U. Then there is a 
constant N such that for any (yi, . . . , y„) G S, if we regard F{yi, . . . , y„, Xn+i) as 
a function of Xn+i on the corresponding U , then it has at most N roots in U if it 
is not identically zero. 

Proof. Since F{xi, . . . ,Xn+i) is algebraic, it satisfies an equation 

d 

(4.1) ^pi(a;i,...,x„+i)X^ = 0, 

with po, ■ • • ,Pd polynomials (not all identically zero). Let N — max^degpi. We 
will show that N has the desired property, using induction on n. 

We can choose the coefhcients pi so that they have no (non-constant) common 
factor. Fix yi e C, and let S" = {(y2, ...,?/„)€ C"-i : {yi, . . .,y„) € S}. Define 
G{x2, . . .,Xn+i) = F{yi ) . Since the coefficients were taken to have no 

common factor, they do not all vanish when we set Xi = yi. Their degrees do not 
increase when we set Xi = yi (or when we remove common factors), so our lemma 
follows by induction on n (applied to G and S"), assuming we can prove it in the 
case n = 0. 

Suppose n = 0. Assuming F is not identically zero, we can divide (13) by 
some power of X to get an equation satisfied by F with non-zero constant term, 
say ph{xi). (A priori, F will satisfy the new equation only where F is non-zero. 
However, since F is holomorphic on [/, its zeros are isolated. By continuity, it 
satisfies the equation at its zeros as well as elsewhere.) Then any root of is a 
root of ph, so F has at most degp/i roots, and hence at most N roots. □ 



of the sum (3.1 
(PropositionllC 



Lemma 15. The exponential sums 

exp(2i <&(£, m- k;n- k)) 

k£l 

remain bounded (uniformly in I) as n goes to infinity, where I can be any subinterval 
of Se, as long as \£\ + \m\ < (1 — 6)n for some fixed S > 0, and e is small enough 
compared to S. 

Proof. To prove this, we will apply the Kusmin-Landau Theorem (Theorem p[). 
Lemma |^ says that d^^{£, m — k; n — k)/dk'^ satisfies the conditions of Lemma ^ 
so there is an absolute upper bound for the number of roots that it can have as a 
function of k while n, and m are held fixed (unless it is identically zero for those 
values of n, and m). Before we apply the Kusmin-Landau Theorem, we break 
Se up into a bounded number of subintervals on which (i$(£, m — k; n — k)/dk is 
monotonic. 
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We have to look at the behavior of d^{£,m — k;n — k)/dk. As in Lemma ^ 
set k = nn, £ = xn, and m — yn. Lemma ^ says that as n goes to infinity, 
d^{£, ni — k;n — k)/dk equals 



(4.2) 



6[~x + y-K, v/(l - k)2 - 2(x2 + {y-nf)^ - 
(l - K - 2x, v/(l - - 2 (x2 + (y - k)2)) +o{^ 



We would like to show that when divided by tt, (E^) stays away from integers. 



After (4.2) is divided by tt, the only possible integral values it can take on are 0, 
±1, and ±2 (assuming n is large enough). If we ignore the 0{\/n) term, the rest of 
the formula is the difference of the arguments of two points on the same horizontal 
line (divided by tt). Thus, it cannot be ±2. It can be only if the points coincide or 
are on the horizontal axis. It can be ±1 only if the points are on the horizontal axis. 
The points coincide x + y = 1, which is impossible (since \t\ + |r7i| < (1 — 5)n). 
They are on the horizontal axis iff 

(4.3) a;2 + (y - k)2 ^ (1 - k)2/2. 



The definition (2.8) of implies that 

x'' + {y-Kf < (l-£)(l-«)2/2. 



so no fc € S'e gives a k satisfying ( [4.3[ ). (Note that k = 1 is impossible since then 
|x| + M = |0| + |l|>l-<5.) 

In fact, the above argument, combined with continuity considerations, shows 
that the two points cannot get arbitrarily close to each other or the horizontal axis, 
and they clearly cannot get arbitrarily far from the origin. Thus, even taking into 
account the 0{l/n) term, {d^{£, m — k;n — k) / dk) / tt really does stay slightly away 
from integers as n — *■ oo. Hence, the Kusmin-Landau Theorem tells us that the 
exponential sums are bounded (uniformly in /). □ 



Proposition 16. The sum 

^ exp(2z <&(£, k;n~ k)) 

^V(n-fc)2-2£2-2(m-fc)2 

is 0{n^^) as n goes to infinity, as long as \£\ + \m\ < (1 — for some fixed 6 > 0, 
and e is small enough compared to S. 

Proof. Let 

aik) = = 

7rV(rt-fc)2-2£2„2(m-fc)2 

and 

fe-i 



h{k)= ^ exp(2i$(£,m-fc';n-fc'))- 



k'=k„ 



For k G Se, a{k) ~ 0{n ^) (by Lemma ||) and b{k) is bounded (by Lemma |T5|) 
Suppose \b{k)\ < B for aU k e Se- 
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To bound the sum in the statement of the proposition, we wiU apply summation 
by parts. We have 

exp(2.y,^ fe. fc)) ^ 

fees, '^V (" - - _ 2(m - ^.^^^ 

= ^ a(fc)6(fc + 1) - ^ a(fc)6(fc) 

a{k~l)b{k)- a{k)b{k) 

fc=fc„i„ + l k=k,„i„ 

^ 6(fc)(a(fc-l)-a(fc)) + 0(n-i). 

fc=fe„i„+i 

This sum is bounded in absolute value by Bj^k ~ 1) ~ "^^C^)! +0{n^^). The 
function a{k) is monotonic on (— cx),2m — n) and (2m — n,oo) (on the subin- 
tervals where it is real, of course), so within each of these intervals, the sum 
\a{k — 1) — a{k)\ telescopes. The boundary terms are 0{n~^), and hence the 
entire sum is 0{n^^). □ 

5. Conclusion of the Proof 

The results proved in the preceding three sections give us Theorem |2|, a weak- 
ened version of Theorem |l|, in which we are restricted to estimating the placement 
probabilities at normalized locations {x,y) with |a;| + |?/| < 1 — 5 for some fixed 
S > 0. That is, we are required to keep {x,y) from getting too close to the bound- 
ary of the diamond. Here we will show how the restriction on [x, y) can be relaxed, 
provided that we are careful to stay away from the points (±^, \ )- 

Fix (5 > 0, and consider the region in the Aztec diamond of order n defined (rela- 
tive to normalized coordinates) by the constraint x^ +y'^ > ^ + (5. This region splits 
up into four pieces. Proposition |l^ tells us that the north-going placement proba- 
bilities tend uniformly to 1 in the northern piece and to in the other three pieces. 
The only regions that are not covered by this method are four small curvilinear 
trapezoids near the points (±i,±i), defined by the inequalities x^ + y"^ < \ + 5 
and 1 — (5 < -I- |?/| < 1. If {x,y) stays away from these four points as n goes to 
infinity, then we can indeed conclude that the placement probabilities for north- 
going dominos at location (x,?/) are as claimed in Theorem £1 This completes the 
proof of the main theorem, except near the points (±i,— 2), which we will now 
deal with. 

Let R be the subregion of the Aztec diamond of order n consisting of the two 
lower of the four curvilinear trapezoids defined hy x^ + y"^ < 5 + <5 and 1 — 6 < 



\y\ < 1. In R, we use the inequality (1.4). It says that for /i > 0, 



F{£, to; n) < m + h;n + h). 

If n is sufficiently large, then for each point (^, to) in i?, there exists an h such 
that the point {i, m + h) of the diamond of order n + h has normalized coordinates 
satisfying x^ + y"^ < ^ + S, 1 - 26 < \x\ + \y\ < 1 - S, and y < 0. Let S be the 



set of all {x,y) satisfying these three constraints. Inequality (1.4) tells us that the 
placement probabilities within R are at most as large as those within S. However, 
the part of Theorem |l| that we have already proved gives estimates for the placement 
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probabilities in S, and shows that they tend uniformly to as (5 —> 0. (To see the 
convergence to most easily, look at the level curves of the placement probabilities.) 
We thus conclude that as i5 ^ 0, the placement probabilities in R tend uniformly 
to 0. This completes the proof of Theorem |^. 

Unfortunately, our techniques do not give us an explicit bound for the difference 
between the placement probabilities and the arctangent formula in an Aztec dia- 
mond of a given order. This is not because the methods are inherently ineffective; 
rather, it is because we have not determined the dependence on e in the 0{n~^) 
term of the error bound in Theorem |l2|. To determine it, we would have to do so 
for the error term in Proposition ^, which seems more trouble than it would be 
worth (but could perhaps be done). 

Using these techniques, we can also prove the arctic circle theorem of | JPS|. One 
direction, that the regions outside the inscribed circle are indeed frozen, follows from 
Proposition ^ To see this, consider the region R defined (relative to normalized 
coordinates) hy + y"^ > + e, with £ > 0. The number of domino spaces in this 
region of an Aztec diamond of order n is less than n^, so the probability that a 
non-north-going domino will appear in the subregion with y > ^, or that a north- 
going domino will appear in the rest of R, is exponentially small, by Proposition ^ 
From this, we see that with probability approaching 1 (as n goes to infinity), all 
the dominos in R will be aligned in brickwork patterns, and thus contained in the 
polar regions. This is half of the arctic circle theorem. 

The other direction, that the polar regions almost never extend substantially 
into the interior of the inscribed circle, requires an additional result for its proof. 
Intuitively, it follows from our main theorem, which tells us that inside the inscribed 
circle all four types of placement probabilities are positive. This trivially implies 
that the polar regions cannot almost always cover a given part of the interior of 
the circle, but showing that they almost never do is harder. We will prove it in 



subsection 6.4 



6. Consequences of the Theorem 

6.1. Height functions. Height functions for domino tilings, which were intro- 
duced in the mathematics literature in ^ (and independently in a slightly different 
form in the physics literature in [Q ) , are a very useful device in the study of tilings 
of simply-connected su bsets o f the plane. (A more general approach to height func- 



tions can be found in |STCR| ].) In any such region that can be tiled by dominos, 
the number of enclosed white squares and the number of enclosed black squares 
under an alternating coloring of the squares must clearly be equal. It follows that 
if one travels around the boundary of the region counterclockwise, then one will see 
a black square on one's left half the time and a white square on one's left half the 
time; to see why, notice that the edges of the square grid that lie within the region 
pair sides of black squares with sides of white squares, so the excess of unpaired (i.e., 
boundary) sides of black squares over unpaired sides of white squares is four times 
the excess of black squares over white squares. As one travels around the boundary, 
any temporary excess of one kind of square over the other that is observed along the 
way represents a "debt" that will eventually have to be paid. Moreover, the same is 
true simultaneously for all the boundaries of all the simply-connected regions that 
are formed by suitable subsets of the tiles in question. Height functions provide a 
uniform framework for keeping track of all these debts simultaneously. 
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Figure 3. A height function for a domino tiling. 



If R denotes a finite, simply-connected region composed of lattice squares that 
have been alternately colored black and white, a height function on R is an integer- 
valued function h on the vertices of the lattice squares which satisfies the following 
two properties for adjacent vertices u and v: first, if the edge from w to u is part of 
the boundary of i?, then \h{u) — h{v)\ = 1, and second, if the edge from it to i; has a 
black square on its left, then h{v) is either h{u) + l or h{u) — 3. It is not hard to show 
that such a function necessarily satisfies a discrete Lipschitz condition: if vertices u 
and V are at distance d from each other in the sup-norm, then \ h{u) — h{v)\ < 2d+l. 
Note also that if h{-) is a height function, then so is h{-) + C for any integer C . 

Every height function on R determines a domino tiling of i?, consisting of domi- 
nos that occupy all the domino spaces that are bisected by edges uv with the 
property that \h(u) — h{v)\ = 3. Conversely, every domino tiling of R arises in 
this way from a height function on R that is unique modulo addition of a global 
constant. We can remove this ambiguity by constraining a particular vertex on 
the boundary of R to have some particular integer value as its height; then every 
domino tiling of R has a unique height function subject to this constraint, and what 
is more, all these height functions agree with one another on the boundary of R. 
For instance, in the case of the Aztec diamond of order n, we set things up so that 
the middle vertex on the west edge of the diamond has height and the middle 
vertex on the northern edge has h eight 2n . (Note that this differs by 1 from the 
height function convention used in [ EKLP| . ) Then the heights for a typical domino 
tiling of an Aztec diamond are as shown in Figure ^. (The shading convention for 
the lattice squares is the same as that in Section]^, i.e., so that the leftmost square 
of each row in the top half of the diamond is white.) 

One can also develop an analogous theory of height functions for other sorts of 
tilings, for example, tilings of regions in the triangular lattice by lozenges (two unit 
equilateral triangles joined along an edge). This theory is simpler geometrically than 
that for domino tilings; for the details, see (See also [BH] for an independent, 
earlier development of height functions for this lattice in physics.) Height functions 
can furthermore be applied to th e sq uare ice model studied by Lieb, as is shown in 
|vB[. The results of subsections 6.2 and 6.3 generalize straightforwardly to other 
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sorts of height functions. However, because the focus of this article is on domino 
tihngs, we wiU not go into the details of the generalization. 

Suppose u, V, and w are three consecutive vertices along a path in a simply- 
connected region R that is tiled by dominos, such that neither the edge uv nor the 
edge ?;w bisects a domino. Then \h{u)^h{v) \ — \h{v)~h{w)\ — 1, with h(w) = h(u) 
if the three points are coUinear and hiw) — h{u) ± 2 if they are not. This principle 
makes it fairly easy to go through the tiling, assigning heights to the vertices. 
Alternatively, one can use this method just to find the heights along the boundary, 
and then find the heights in the interior by the following procedure. To determine 
the height of a particular vertex in the interior of a tiled region, start at the point 
on the boundary of the region due north of the vertex (whose height is independent 
of the tiling) and proceed downward, subject to the following rule: when one travels 
southward along an edge that bisects a north-going (resp. south-going) domino in 
the tiling, the height decreases (resp. increases) by 3, whereas, when one travels 
southward along an edge that bisects a north-going (resp. south-going) domino 
space that is not occupied by a domino in the tiling, the tiling, the height increases 
(resp. decreases) by 1. A similar rule can be formulated for describing how the 
height changes as one travels horizontally through the interior of the diamond. The 
fact that these rules are consistent with each other is a consequence of the fact that 
any region that can be tiled by dominos must contain exactly equal numbers of 
black and white squares. 

If we take the average of all the (finitely many) height functions associated with 
the different tilings of a region, we get a real-valued function on the vertices called 
the average height function. As a consequence of the rule described in the preceding 
paragraph, one can give a simple description of how the average height changes from 
vertex to vertex, in terms of the placement probability p associated with the domino 
space that is bisected by the edge that connects the two vertices. For instance, if 
u and V are neighbors, with u to the north or west of v, then the average height at 
V is equal to the average height at u plus 4p — 1 if edge uv bisects a south-going 
or west-going domino space, while the average height at v is equal to the average 
height at u minus Ap—1 if edge uv bisects a north-going or east-going domino space. 

Here we are interested in the asymptotic behavior of the average height function 
for domino tilings of the Aztec diamond. The arctangent formula tells us that 
these probabilities p are slowly varying, so the average height function is locally 
approximated by functions of the form ax-\-hy + c (with a, 6, and c slowly varying). 
We call the pair (a, b) the tilt of the plane z = ax + by + c. Let us normalize 
our height functions by dividing through by n, both in the domain and in the 
range. Thus, in the limit we expect to see some sort of function 7i(-, •) on {(x,y) : 
< 1} satisfying the piecewise- linear boundary condition H (a;, y) — l—x^+y'^ 
for \x\ + |y| = 1, as well as a Lipschitz condition with constant 2 relative to the 
sup-norm distance. In addition, the formulation in the previous paragraph of how 
the average height changes when moving between vertices tells us that we should 
have ^ = 2{pn-ps) and ^ = 2(p„ -pe), where pn, Ps,Pe, and p^ are the north- 
going, south-going, east-going, and west-going placement probabilities at (a;,j/), 
respectively. 

It can be s hown ( although we do not prove this here) that the domino shuffling 
algorithm of [EKLP] can be interpreted directly in terms of height functions, and 
that half of the values of the average height function for the diamond of order n + 1 
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are equal to certain corresponding values of the average height function for the 
diamond of order n. Hence the average height functions for the Aztec diamonds of 
orders n and n + 1 cannot be too far apart. However, such considerations are not 
sufficient to yield a proof that the normalized average height functions converge to 
a continuum limit. 

The arctangent formula gives us the strength we need in order to conclude that 
a limit exists. Recall that the average height function can be derived by taking 
cumiilativc sums and differences of local placement probabilities, with various co- 
efficients. Taking this assertion to the limit as n ^ oo, and using the fact that the 
placement probabilities approach a continuum limit, we see that the normalized 
average height function also approaches a limit. (It is true that the errors in the 
placement probabilities are going to add, and that there are more and more of them 
as n gets large, but each individual error is going to be small, so that when we divide 
by n the normalized error is small as well.) The limit must be some function 7i(.T, y) 
(defined for |x| + \y\ < 1) with the property that ^ = 2P(y, —x) — 2P{—y, x) and 
^ = 2P{x, y) — 2P{—x, —y) for all {x, y) in the interior of its domain. (Here, 
we have expressed the placement probabilities near (x, y) for all four domino ori- 
entations in terms of V via rotational symmetry.) That is, the tilt of the tangent 
plane at a point (associated with the average height function) can be expressed in 
terms of the local placement probabilities for random domino tilings. This means 
that we ought to be able to reconstruct the function H(-,-) from Theorem 1 via 
integration, making use of the known boundary conditions satisfied by H. If we do 
this, it turns out that H can be written in closed form, and indeed, a formula for 
H can be written especially compactly if one makes use of the formula for P{-,-) 
itself. Specifically, one can verify that the following formula for ■^(•, •) holds: 

Proposition 17. The normalized average height functions for large Aztec dia- 
monds converge uniformly to 



n{x, y) = 2 {yP{x, y) - yP{-x, -y) + (1 - x)P{-y, a;) + (1 + x)P{y, -x)) . 

Proof. We simply check that this formula satisfies the differential equations and 
boundary conditions. □ 



Within the temperate zone, the average height function is real analytic; in each 

of the polar regions, it is an affine function of x and y. On the arctic circle itself, 
away from the points (±5, ±5), the function ?i(a;, y) is differentiable but not twice- 
differentiable. It takes the value 2 at the points (0, ±1) and the value at the 
points (±1,0), with piecewisc-linear behavior on the boundary of the normalized 
diamond. The level set 71 = 1 consists of the two line segments joining midpoints 
of opposite sides of the normalized diamond. 

One may ask, for diamonds of finite order n, how closely the distribution on 
height functions is clustered around its mean value. We will see in the next sub- 
section that the standard deviation of the height at any fixed location in the Aztec 
diamond of order n is at most S^/n. However, numerical evidence suggests that, at 
the center of the diamond, the standard deviation of the height is much smaller — 
more like logn, or perhaps even less than that. 
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Our formulas for ^ and ^ , in combination with the arctangent formula, yield 
(within the temperate zone) the equation 

d'^n d^n 8 



dif dx^ ttVI- 2x2 -22/2' 

We can rewrite this equation in a slightly more illuminating way. For t > and 
+ \y\ < t, define 

H{x,y,t)^m{x/t,y/t). 
That is, we undo the scaling introduced in Section 0. Then off the arctic circle we 
have 

8u[x,y,t}, 



where 



dy^ dx^ 

if a;2 + y2 < ^2/2, and 



u(x,y,t) = <^ ^^/t^-2x^-2v^ 

[O if x2 +y2 > ^2/2. 

This function is a fundamental solution to the wave equation in two dimensions, 
with speed of propagation 1/V2- That is, u is a distribution satisfying 

u{x,y,0) 0, 

and ^ 

— {x,y,0) = S{x,y), 

where S is the (two-dimensional) Dirac delta function. (For more details on funda- 
mental solutions to the wave equation, see |^ p. 164].) 

Note that Proposition ^ shows that, except for an oscillating factor, the creation 
rates also behave like 2u. William Jockusch has shown in personal communication 



how to use a generating function developed in ||GIF(] to explain this behavior, by 
viewing the creation rates as numerical approximations to a solution of the wave 
equation. 

He has also pointed out that from his methods, one ought to be able to deduce 
a weak version of Theorem ^ More specifically, one should be able to show that 
in any macroscopic subregion of a randomly tiled Aztec diamond of order n (i.e., 
any subregion of size comparable to that of the diamond) , the expected density of 
north-going dominos is within o(l) of that predicted by integrating the arctangent 
formula; in particular, this would suffice to prove Proposition Unfortunately, 
his methods would not rule out the possibility of small-scale fluctuations in the 
placement probabilities, such as one gets if one looks at placement probabilities for 
all horizontal domino spaces rather than just the north-going ones. 

6.2. Robustness. The formula for the average height function that was derived 
in the preceding section from the arctangent formula applies not only to Aztec 
diamonds, but also to all regions that approximate them in a suitable sense. (Here, 
as hereafter, the term "region," without qualifiers, should be understood to refer to 
finite regions in the plane that are unions of lattice squares and can be tiled with 
dominos.) It is not enough that the region being tiled should have a boundary that 
is roughly "Aztec" in shape. For instance. Figure ^ shows a random domino tiling 
of a region obtained from the Aztec diamond of order 32 by adding some squares 
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Figure 4. A modified Aztec diamond with far more tilings. 



along its boundary, while the region that is shown in Figure ||, also studied in |SZ|, 
was obtained by adding an extra row of length 64 in the middle of the diamond. 
(These random tilings were obtained via the method described in [PW| and are 
indeed truly random, to the extent that pseudo-random number generators can be 
trusted.) In neither case do we get behavior consistent with the arctangent formula. 
On the other hand, Figure ^ shows a random tiling of an Aztec diamond to which 
two rows of length 64 have been added, and the resemblance to Figure]^ is evident. 

The sense of mystery dissolves if one considers the behavior of the height function 
along the boundary in each of the three cases. In the first case, the height function 
is nearly constant along the boundary; in the second, the direction of change of 
the height function is the same along the southwest and northwest edges (and also 
the same along the southeast and northeast edges); and in the third, the direction 
of change of the height function switches as one rounds any of the four corners of 
the region. Since it is the third situation that resembles the boundary behavior 
of height functions of Aztec diamonds, it is not surprising that the third situation 
should also give behavior in the interior that is similar to what one sees for Aztec 
diamonds. (In fact, if we view the third region as an Aztec diamond of order 33 with 
two vertical dominos removed, then since almost all tilings of the Aztec diamond 
contain those two dominos, the local statistics in the third region differ little from 
those in the Aztec diamond of order 33.) 

Note, incidentally, that if instead of adding a row of length 64, as we did in 
Figure ^, we removed a row of length 64, then the resulting region is easily seen 
to have a unique tiling, consisting entirely of horizontal brickwork. Although this 
situation may seem trivial, it can shed some light on what is happening in Figure |[ 



28 



HENRY COHN, NOAM ELKIES, AND JAMES PROPP 



Figure 5. A modified Aztec diamond with far fewer tilings. 



The horizontal brickwork pattern seen almost everywhere in Figure ^ is the unique 
arrangement of dominos such that the height increases (or decreases, depending 
on whether the dominos in the pattern are north-going or south-going) as quickly 
as possible as one moves vertically. In an Aztec diamond with a row of length 64 
removed, the heights on the boundary are such that the only way to span the gap 
between the heights on the lower half of the boundary and those on the upper half 
is to change at this rate. Thus, the only way to tile the region is with a brickwork 
pattern. In the case of Figure ||, the occurrence of an extra row of length 64 gives 
the height function a tiny bit of slack, and we can see where this slack gets used by 
following the fault-line that runs from left to right. 

Of course, we could have predicted ahead of time that small modifications of the 
shape of the boundary can have a drastic impact on the tiling situation, since for 
instance adding a single square to a region (or removing a single square) can create 
a region with odd area, which cannot be tiled at all. Hence, we will want to assume 
that all the regions we discuss actually admit tilings by dominos, as stipulated in 
the first paragraph of this subsection. 

We will show in this subsection that regions similar to Aztec diamonds, such 
as Figure ^, have approximately the same average height functions as the Aztec 
diamonds they resemble. This will follow as a consequence of a more general result, 
asserting that the value of the average height function depends in a continuous 
manner on the values of the fixed heights along the boundary. That is, if one 
modifies the shape of the boundary in such a way that the height function along 
the new boundary, when plotted in three dimensions (the two original dimensions 
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Figure 6. A modified Aztec diamond with Aztec-like height function. 

plus a third dimension for height), is close to the graph of the old, the average 
heights of vertices in the interior will not change very much. 

Before we can do this, we first prove a general monotonicity result about height 
functions. The idea for this approach was suggested by Robin Pemantle in personal 
communication. Let 7? denote a simply-connected region in the plane with some 
fixed checkerboard coloring, and let V be the set of vertices in R. Let V be a 
subset of V that contains all the vertices on the boundary of R; we will assume 
that V is connected, in the sense that the subgraph of the square grid induced by 
the vertex set V is connected. A partial height function is a function f : V ^ Z 
subject to the local constraint that if u and v are adjacent vertices such that the 
directed edge from w to w has a black square on its left, then f{v) — f{u) is either 1 
or —3. It is called complete height function if it is defined on all of F; a complete 
height function / extends a partial height function / if it agrees with / where / is 
defined. 

Throughout this subsection (and the next), H will denote a complete height 
function chosen at random (according to some distribution); thus, for any vertex 
V, Exp[H{v)] (the expected value of H{v)) is the value at v of the average height 
function. 

Given a connected subset V of V that contains all the boundary vertices, and 
given a partial height function / on V', we let /i/ denote the uniform distribution 

on the sot of complete height functions that extend / to V. 

Lemma 18. /// and g are two partial height functions defined on V and agreeing 
modulo 4, mth f < g, then is stochastically dominated by ^g. That is, there 
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exists a probability measure tt on the set of pairs {f ,g) of complete height functions 
extending f and g respectively, such that 

Y,Afrg) = Hifl 

g 

/ 

and 

t({(/,.9):/<.9}) = 1. 

Proof. We use induction on the size of V \V' (holding V fixed and varying V'). 
The case where this set is empty is triviaL Assume that the lemma is true whenever 
\V \ V'\ = k — 1, and suppose we have a situation in which |y \ = fc. It clearly 
suffices to consider the case in which f{v) < g(v) for some vertex v in V' that is 
adjacent to at least one vertex in V\V' . Let w be a vertex in W = V\V' adjacent 
to V. 

Given that f{v) has some specific value, any extension of / to V" — V U {w} 
would have to give w height /i or ft, — 4 (for some particular h whose value we don't 
care about — it's f{v) plus or minus 1 or 3), while any extension of g to V" would 
have to give w height h' or h' — 4 (with h' determined from g(v) the same way h 
is determined from f{v)). Because / and g agree modulo 4 on V' and h' > h, we 
have h' ~ A> h. 

Let /{ and be the two extensions of / to V" that assign w height h and /i — 4, 
respectively, and let g[ and g2 be the two extensions of g to V" that assign w height 
h' and h' — 4, respectively. (If such extensions do not exist, it is not a problem, 
as we will see below.) The distribution /x/ is a weighted superposition of fj,f^ and 
fifi^, where the ith term (i = 1 or 2) is given weight proportional to the number 
of extensions of f^ to V (which should be taken to be zero in the case where the 
extension to V" does not exist). Similarly, Hg is a superposition of /igj and /ig^. 
Since < g'j for all i,j in {1,2}, and h = h' (mod 4), we can use our induction 
hypothesis to conclude that ^.f is stochastically dominated by ^g', for all i,j, which 
implies that /i/ is stochastically dominated by /Xg, as was to be shown. □ 

Corollary 19. // / and g are two partial height functions on R defined on V' and 
agreeing modulo A, with f < g + A, then for all v, Exp[i/(i;)] under the measure /z/ 
is at most 4 more than Exp[i7(w)] under the measure fig. 

Proof. Apply Lemma |l^ to the partial height functions / and 5 + 4. □ 

For applications of Lemma |l8| and Corollary it is important to note that the 
values of height functions on connected regions are determined modulo 4, given the 
value at any one point, because the defining conditions for a height function imply 
that if two height functions agree modulo 4 at any point, then they do so at each 
neighboring point. Also, given two partial height functions defined on different 
sets, we say that they agree modulo 4 if all height functions extending them agree 
modulo 4. 

Proposition 20. Suppose that Ri, i?2 are two simply- connected regions in the 
plane, with mandated partial height functions fi, f2 along their boundaries that 
agree modulo 4, such that every vertex vi on the boundary of Ri is within sup- 
norm distance Ai of some vertex V2 on the boundary of R2, and vice versa, and 
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such that whenever vertices Vi and V2 on the respective boundaries are within sup- 
norm distance Ai of each other, the heights ,fi{vi) and f2{v2) are within A2 of each 
other. Then, for any v in RiH R2, the expected value of H(v) under /i/j^ and the 
expected value of H(v) under differ by at most 2Ai + A2 + 1. 

Proof. Let /i.max be the highest extension of /i to Ri, let /i.min be the lowest 
extension of /i to and define /2,max and /2,min similarly. (It is not hard to show 
that the complete height functions extending a given partial height function form a 
lattice under the usual partial ordering, so it makes sense to talk about the highest 
and lowest extensions.) Let v be on the boundary of Ri fl R2 (and hence on the 
boundary of Ri ot R2). If u is on the boundary of then we can find a nearby 
w on the boundary of R2 so that 

/l,max(w) = fl{v) 

< /2(W)+A2 

= /2,min(w) + A2 

< /2,mi„(v)+2Ai + l + A2, 

while if V is on the boundary of R2, then wc can find a nearby w on the boundary 
of i?i so that 

/l,max(w) < fl ,inax 

(w) + 2Ai + 1 
= /i(w) + 2Ai + l 

< /2(i;) + A2 + 2Ai + l 

= /2,mi„(w)+A2 + 2Ai + l. 

Since the two height functions agree modulo 4 at v, /i,max(i^) < /2,min(i') + 'iK, 
where is the greatest multiple of 4 that is less than or equal to 2Ai + A2 + 1. 
It follows from this (and the corresponding inequality in the other direction) that 
if /{ is any extension of /i to i?i and /2 any extension of /2 to R2, then for each v 
on the boundary of i?i fl i?2, fi{v) differs from /2(w) by at most AK. 

Now let V be any vertex in i?i fl R2. If we compute Exp[_ff (w)] by conditioning 
on the heights on the boundary of RiD R2, then it follows from Corollary ^ that 
the expected value of H{v) under /i/^ differs by at most AK (and hence at most 
2Ai + A2 + 1) from its expected value under /i/^. □ 

As an application of this result, we may consider a modification of the Aztec 
diamond of order n, whose symmetric difference with the true Aztec diamond of 
order n is a narrow fringe around the border of the true diamond, of width o{n). 
Suppose that the black and white squares of the symmetric difference are equinu- 
merous, and moreover that they are not segregated but intermix in such a manner 
that the height function along the border of the modified diamond is within o(n) 
of the height function along the border of the true diamond. Lastly, let us suppose 
that the modified diamond has at least one domino tiling. Then we can conclude 
that the average height function for the modified diamond is within o{n) of the 
average height function for the true diamond. 

Notice that these results give us no direct information about how individual 
placement probabilities change in response to small changes in the shape of the 
boundary, though some weak information can be obtained by way of the height 
function. It would be quite interesting to obtain robustness results for the placement 
probabilities themselves. 
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6.3. Variance. In Proposition 17 of subsection 6.1 we gave a formula for the nor- 



malized average height function, or rather its limit as the size n of the Aztec di- 



amond goes to infinity, and in subsection 3.2 we showed that the same formula 
applies to many regions that are roughly similar to the Aztec diamond. However, 
we obtained no information about how closely a typical height function for a region 
(an Aztec diamond or something else) should approximate the average height func- 



tion. Here we use Azuma's Inequality | AS, p. 85] to bound the amount of variation 



that values of random height functions are likely to exhibit. 

Let H denote the (unnormalized) height function corresponding to a random 
domino tiling of some simply-connected region in the plane, and let v denote a 
vertex in the region, such that there is a path of m vertices from the boundary of 
the region to v. We will show in this subsection that the variance of the random 
variable H{v) is at most 64m. In fact, we actually get a stronger result: 

Theorem 21. Let f be a partial height function defined on the boundary of a 
simply- connected region R, and let v be a vertex in the interior of R, such that 
there is a path of m vertices from the boundary of R to v. Then, for all c > 0, 
the probability that H(v) (the value of a random height function at v under the 
uniform distribution fif) differs from its expected value by more than C'^/m is less 
than 2e-'='/32. 

Proof. Let xo,xi, . . . ,Xm-i = f be a lattice-path connecting a point Xq on the 
boundary of R to the point v. Let Fk be the partition of the space of possible 
height functions in which two height functions are regarded as equivalent if they 
agree at xo, xi, . . . , Xk-i- Let Mk be the conditional expectation Exp [iJ(ti the 
function from the set of height functions to the reals that assigns to each height 
function h the average value of h'{v) as h' ranges over all height functions in the 
equivalence class of h. 

Note that Mm is just the function H{v) itself, while Mq is the average value of 
the height at v, averaged over all height functions. The fimctions Mq, Mi, . . . , Mm 
form a martingale; that is, 

Exp[Mfc+i|Ffe] 

On each component of Fk, Mk+i — Fjxp[H {v)\Fk+i] takes on at most two distinct 
values, according to the two different values of H(xk) that are consistent with the 
already-known values of H{xq),H{xi), . . . , H{xk-i)- From Corollary |l^, we see 
that these two values of Mk+i differ by at most 4. Since Mk is their weighted 
average, it follows that Mk and Mk+i never differ by more than 4. Then, applying 
Azuma's Inequality (Corollary 2.2 on page 85 of | |AS| |) to the quantities Mfc/4, we 
get 

Prob[|M„ - Mol/4 > ty/7^] < 26"*'/^ 
Replacing t by c/4, we get 

PToh[\H{v) - Exp[i/(D)]| > cVm] < 26^^='/^^. 
This completes the proof. □ 



If we are interested in estimating the variance, we can derive a consequence of 
the preceding inequality: Assuming for simplicity of derivation (and without loss 
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of generality) that Exp[77(u)] — 0, we have 

Vai[H{v)] = Exp[{H{vyf] 











< 

= 64m 



Prob[(iJ(w))^ > x] dx 
Vioh[\H{v)\ > ^/x\dx 

2g-x/(32m)^^ 







As a final aside, we mention that our proof of Theorem |2l| also yields the following 
more general result: 

Proposition 22. Let R be a simply- connected region in the plane and let v,w be 
vertices in the interior of R, such that there is a path of m vertices from v to w, 
staying entirely within R. Then, for all c > 0, the probability that H{v) — H{w) 
( under the uniform distribution on domino tilings of R) differs from its expected 
value by more than C\pm is less than 2e^^ Z"^^. 

6.4. The arctic circle theorem. We will now use Theorem pi| to complete the 
proof of the arctic circle theorem, which we began in Section 0!^We still need to 
show that the polar regions almost never extend very far into the interior of the 
inscribed circle +?/^ = ^ (defined relative to normalized coordinates). Let e > 0, 
and consider the region R in an Aztec diamond of order n defined by a;^+y^ < \ — e. 

A domino is in the north or south polar region if and only if the heights on the 
vertices of the domino are equal to those at the same locations in the all-horizontal 
tiling, which is the minimal tiling of the Aztec diamond (under the partial order- 
ing of tilings induced by comparison of height functions). An analogous statement 
connects the other two polar regions to the all-vertical tiling, which is the maximal 
tiling. Proposition |l^ shows that, asymptotically, the average height function dis- 
agrees with the minimal and maximal height functions within the inscribed circle 
(although outside of that circle it agrees with one or the other). In particular, 
if a domino in R is part of the polar regions, then the heights on it differ from 
the average heights at those locations by an amount at least proportional to n. 
(Of course, the constant of proportionality depends on e.) We see that, by taking 
c = \/n in Theorem |2^, the probability that a domino in R will be part of the polar 
regions is exponentially small in {^/n)^ = n. Since the number of dominos in R is 
on the order of n^, the probability that any will be contained in the polar regions 
is exponentially small. 

We have now proved a slightly stronger version of the arctic circle theorem than 



that proved in [JPS|. There, it is shown that for any e > 0, for sufficiently large 
n, the boundary of the temperate zone stays within distance en of the arctic circle 
with probability greater than 1 — s. We have shown that this probability differs 
from 1 by an amount exponentially small in n. 

6.5. Heterogeneity. The arctangent formula gives us an indication of a certain 
sort of local homogeneity: places in the tiling that are close together tend to be 
governed by the same statistics. Here we would like to prove a converse result, and 
show that within the temperate zone, places in the tiling that are far apart tend 
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to be governed by different statistics. More precisely, we would like to show that 
within the temperate zone, the quadruple 

{Pn,Pe,Ps,Pw) = {'P{x,y),V{-y,x),V{-x,-y),V{y,-x)) 

(whose components are the four placement probabilities near the location {x,y)) 
uniquely determines x and y. This means, to put it somewhat fancifully, that if 
you found yourself stranded somewhere in the temperate zone of a random domino 
tiling of a huge checkerboard colored Aztec diamond, then, provided that you had 
a compass to tell you which way was north, you could determine your relative 
position within the diamond merely by examining the local statistics of the tiling. 

The heterogeneity claim is not hard to prove, since we know that the level sets for 
all four placement probabilities are arcs of ellipses having a very specific geometry. 
In particular, level sets for p„ and Ps are arcs of ellipses that intersect in at most 
two points, and these two points have the same y-coordinate; similarly, level sets 
for pe and Pw are arcs of ellipses that intersect in at most two points, and these 
two points have the same x-coordinate. It follows that if two points have the same 
probability quadruples, they must have the same x- and y-coordinates; that is, the 
two points must coincide. 

However, we wish to prove more. Consider that the elements of the quadruple 
sum to 1, so that the quadruple has three degrees of freedom. However, x and y to- 
gether embody only two degrees of freedom, so as x and y sweep through their range 
of joint allowed values, the quadruple determined by x and y will not sweep through 
the full set of probability vectors of length 4. On the other hand, the asymptotic 



normalized average height function Tl introduced in subsection 6T manifests exactly 



two of the degrees of freedom of {pmPe,PsTPw) in its first-order derivatives. What 
we would like to show is that inside the temperate zone, the map {x, y) ^ ^) 
is one-to-one, and has as its range the region {(s,t) : \s\ -f |t| < 2}. (The possible 
significance of this fact will be explained more fully in Section ^) Putting it dif- 
ferently, we may say that if one views the graph of the restriction of the function 
Ti to the interior of the temperate zone as a surface, then the Gauss map from the 
surface to the sphere is injective. 



To prove the claim, we first note that, as discussed in subsection 3.1, ^ 



2pw — 2pe and ^ = 2p„ — 2ps. With y fixed and x increasing, pe increases while 
Pw decreases, achieving equality (by symmetry) at x = 0. Thus, the sign of ^ tells 
us the sign of x, and similarly, the sign of ^ tells us the sign of y. Hence, to prove 
the injectivity of the map, it suffices to focus on the part of the temperate zone 
that lies in the interior of one particular quadrant, say the second. Within that 
quarter-disk, ^ and ^ are both non-negative functions, taking the values on 
the respective axes a; = and y = and increasing as one moves away from these 
axes. These monotonicity properties do not of themselves rule out the possibility 
that ^ and ^ have the same value for two different points in that quarter-disk, 
so we must resort to a slightly more arduous approach. 
Using the arctangent formula, one can check that 



(? 



W j 1 - a;2 - 3y2 



cos 



( TV d-H 
\2 dx 



1 — Sx^ — y"^ 
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and 

X 



sm 



( TT dH \ 
\2 dx ) 



(If ix^ + = 1, then the first ratio is not defined. However, since x^ + < ^, 
either the first ratio or its reciprocal is defined.) Given the values of these two 
ratios, there are in general at most two possibilities for {x,y), only one of which 
will be in the desired quadrant. The only case in which knowledge of the two ratios 
does not restrict us to at most two possibilities for {x, y) is when the first ratio is 1. 
This happens iff ^ = i.e., along the line through the origin that bisects the 
quadrant. Since one can check using the explicit formulas for ^ and ^ that the 
partial derivatives increase as one moves away from the axes along that line, they 
still determine {x,y). It follows that the map {x,y) i-^ ^) is injective on the 
quadrant, as was to be shown. 

Now we will see that the map is in fact a surjection to the set {{s, t) : \s\ + \t\ < 2}. 
If one sets x = {1 — t — ct^)/2 and y = {1 -\- 1 — ct^)/2 with c > ^ (so that {x,y) 
lies on a parabola that is symmetric about the axis x = y and that lies inside the 
closed temperate zone in the vicinity of (i, i)), then, sending t to zero from above, 
we find that the north-going and east-going probabilities tend towards 

11 1 1 

- + -tan"-^ . 

2 TT V2c- 1 

and 

11 _i -1 
' — tan 



2 TT V2^^ 

which sum to 1 for all c between ^ and infinity and which vary (as an ordered 
pair) over the open segment connecting (1, 0) to (i, i), as c goes from i to infinity. 
Plugging the limits p„ — s- p, Ps — s- 0, Pe — > ^ — P, Pw ^ into the formulas 
^ = 2pw — 2pe and ^ ~ 2p„ — 2ps, we find that the boundary of the open 

square {{s,t) : \s\ + \t\ < 2} consists of limit points of the set of tilts (^, ^) 
that are achieved by the average height function in the temperate zone, and hence 
(by continuity) that that we do indeed obtain the open square as the set of tilts 
achieved by the average height function in the temperate zone. 

6.6. Entropy. The entropy of a random variable that takes on any of N values 
with respective probabilities gi , . . . , gAr is defined as J2iLi ~1i log li (with log = 
by convention); for example, the entropy of a uniform random domino tiling of the 
Aztec diamond of order n is log 2, because there are exactly 2"("+^)/^ tilings 



(see [ EKLP| for a proof). We have seen that for large n, nearly all of this entropy 



is due to the variety exhibited inside, as opposed to outside, the arctic circle. It 
would be good to have more quantitative information on this. Specifically, given 
a patch of an Aztec diamond, one can define a random variable whose values are 
the near-tilings of the patch that result from restricting a uniform random tiling of 
the Aztec diamond to just the patch (such near-tilings are allowed to have untiled 
squares along the boundary of the patch), and one can consider the entropy of this 
new random variable. If the patch is very large (while not long and skinny, for 
example like a 2 x n rectangle), but the order of the Aztec diamond is much larger 
still, then we believe that this entropy, when divided by the area of the patch, is 
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close to a value which we would call the local entropy, and which would depend 
only on the normalized location of the patch. 

In this subsection, we make a small start towards calculating local entropy by 
showing that it vanishes outside the arctic circle and that it is positive inside the 
arctic circle (assuming it is well-defined there). Assuming that local entropy is 
well-defined everywhere, this gives us another way of interpreting the arctic circle, 
namely as the boundary between the zero-entropy region and the positive-entropy 
region. 

The vanishing (and perforce the well-definedness) of local entropy outside the 
temperate zone is a simple consequence of the arctic circle theorem. To prove the 
other half of our claim, consider an m x m patch P sitting inside the temperate 
zone of an extremely large Aztec diamond, with m even. 

If a, 6, c, and d are the northwest, northeast, southwest, and southeast squares in 
a 2 X 2 block in a plane region R that can be tiled by dominos, then the proportion 
of tilings of R that have a horizontal domino covering squares a and b and another 
horizontal domino covering squares c and d (write this proportion as Pab,cd for short) 
is clearly equal to the proportion Pac,bd of tilings that contain vertical dominos 
covering squares a and c and squares b and d; moreover, by one of the lemmas 
proved in ||GIP |, both proportions are equal to PabPcd + PacPbd, where Pab denotes 
the proportion of tilings that have a domino covering a and b, etc. (that is, Pab, Pcd, 
Pac, and pbd are just placement probabilities under uniform random tiling). In our 
particular situation, if one looks inside the patch P taken from the temperate zone 
of a large Aztec diamond, all four placement probabilities are bounded away from 
zero, say by e > 0, so the probability that a random tiling contains a 2 x 2 block 
centered at any particular vertex in P is at least 4e^. In particular, we can look 
at the (m/2)^ vertices that are at the centers of the (m/2)^ non-overlapping 2x2 
blocks into which P can be naturally decomposed. Using linearity of expectation, 
we can see that the expected number of such 2x2 blocks in a random tiling of the 
Aztec diamond is at least rn^e^ . However, this allows us to set a lower bound on 
the entropy, as measured by the variety of configurations one sees locally. For, by 
freely rotating these blocks (i.e., changing horizontal blocks to vertical blocks or 
vice versa), we can create 2™ ^ other local patterns, all equally likely. Standard 
techniques in information theory permit one to conclude that the entropy of the 
near-tiling of P is at least log 2 times the area of P. 



7. Further Results 

Although we have phrased our results in terms of domino tilings, there is an 
easy equivalence between domino tilings of finite regions and dimer configurations 
on certain finite graphs. Specifically, if we replace each square cell by a vertex, and 
draw an edge connecting any two vertices whose associated cells are adjacent, then 
a domino tiling of a region corresponds to a dimer-cover of the derived graph, that 
is, to a set of edges of the derived graph with the property that every vertex belongs 
to exactly one of the chosen edges. In this way, the study of domino tilings is seen 
to be equivalent to the study of dimer-covers, which is one of the better-understood 
statistical mechanics models in two dimensions. Studying domino tilings of special 
regions, such as Aztec diamonds, is tantamount to studying the dimer model in 
the presence of special boundary conditions. The uniformity of the distribution 
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corresponds to a degenerate situation in which all dinier configurations have the 
same energy. 

There has been surprisingly little work on the behavior of the dimer model in the 
presence of general boundary conditions; researchers in statistical mechanics have 
tended to study either toroidal (i.e., periodic) boundary conditions or boundary 
conditions that correspond to domino tilings of a rectangle. Our work can in a 
sense be regarded as a somewhat strange chapter in the study of the dimer model, 
in which highly unphysical boundary conditions are imposed. (Precursors of this 
research include Q], [^], and |S^. 



In his original article on the dinier model [Ka], Kasteleyn considered imposing an 



energy function that favors one orientation of dimer over another (horizontal versus 



vertical). The authors of |GIP] followed this lead, and showed how their methods 
also led to exact results for random domino tilings of the Aztec diamond when the 
distribution was skewed towards dominos of a particular orientation. Here, we will 
state the results that follow from applying the methods of this paper to the case of 
biased tihngs. 

Let p be strictly between and 1. For each n, there is a unique probability 
distribution on the tilings of the Aztec diamond of order n (in fact, on any simply- 
connected region) such that given any tiling of all of the diamond except for a 2 x 2 
block, the conditional probability that the 2x2 block will contain two horizontal 



dominos is p. For more details on this distribution, see [|JPS[| or |GIP[ . We call this 
the Gibbs distribution with bias p. (For more information on Gibbs distributions 
in general, see Q.) 

The main difference between the biased distribution and the uniform distribution 
is the shape of the temperate zone. We will see shortly that, in the biased case, its 

2 2 

boundary is given by the "arctic ellipse" — + jz- ~ 1 (in normalized coordinates). 



It was conjectured in [IPS| that the analogue of the arctic circle theorem holds 
in the biased case. Our methods prove that conjecture, as well as an arctangent 
formula that describes the behavior within the temperate zone. 

We begin by defining the biased placement probabilities Pp(^, m;n) the same 
way we defined the ordinary placement probabilities (except, of course, that we use 
the biased distribution). The biased creation rates are also defined analogously to 
the ordinary creation rates, by 

Cr„(^, m; n) = -(Pp(^, m; n) — Pp(^, m — 1; n — 1)). 
P 



The proofs depend on a biased version of Proposition ^, which is proved in JGIP | 



To state it, we will need a more general form of Krawtchouk polyn omia l. Define 
Cp{a,b;n) to be the coefficient of .z" in (l + (l-p)z/p)"~''(l-z)^ (See [@ p. 151].) 



We need the following result from | GIP 



Proposition 23. Let < p < 1, and set a = {£ + ni + n)/2 and b= {£ — m + n)/2. 
If a,b Z, then 

Cip{£, m;n + 1) = Cp{a, b; n)cp{b, a; n)p". 
Otherwise, CTp(£, m; n + 1) = 0. 

Using this proposition, straightforward modifications to the proof of Proposi- 
tion H prove the following generalization: 
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Proposition 24. Fix e > 0. Ifj + fz^<{l- e)n'^ and i + 711 = 71 (mod 2), th. 
Crp(^,m;n+ 1) = p\ ^ ^ ) + 0^{7i ^) 

TXyJ {p — p'^)7l^ — (1 — p)P — 

for some Junction $p(^, m;n), which can he deternriined explicitly. 

Every result needed for the proof of Theorem ^ (such as the creation rate esti- 
mates outside the arctic ehipse) has a straightforward generahzation to the biased 
case; in the interest of saving space, we wih omit their statements. The proofs are 
completely analogous to the proofs for the uniform distribution. One arrives at the 
following biased counterpart to Theorem 

Theorem 25. Let < p < 1, and let U be an open set containing the points 
(±p, 1 — p). If {x,y) is the nornrialized location of a north-going domino space in 
the Aztec diamond of order n, and [x, y) ^ U , then, as n 00, the placement 
probability at {x,y) for the Gibbs distribution with bias p is within o(l) ofVp{x,y), 
where 

Pp(a;,y) = <jl ^ if^ + ^>landy>l-p,and 



, if2^ + J^<l. 

The 0(1) error bound is uniform in {x,y) (for {x,y) ^U). 

Similarly, the south-going, east-going, and west-going placement probabilities 
near {x,y) are approximated by Vp{—x,—y), Vi-p{—y,x), and Vi~p{y, —x), re- 
spectively. This follows from Theorem ^ by rotational symmetry. 

One can also prove biased versions of the robustness and variance results from 



subsections 6.2 and 6.3. (In fact, the proofs are practically identical to the proofs 



given in those subsections.) Using them in combination with the same methods 



used in subsection 6.4, we can prove a slightly strengthened version of the "arctic 



ellipse conjecture" from |JP5] 



Theorem 26. Let < p < 1, and e > 0. The probability that, in a random domino 

tiling with bias p of an Aztec diamond of order n, the boundary of the polar regions 

2 2 

is more than a distance e in normalized coordinates from the ellipse ^ -I- = 1 
is exponentially small in 71. 

8. Speculations 

In this article we have focused primarily on one particular family of finite regions, 
namely, Aztec diamonds. Here we will indicate what it might mean to have a theory 
that would apply to all simply-connected finite regions, and how Aztec diamonds 
might play a role in the project of classifying the different possible local behaviors 
that random tilings of such regions can exhibit away from their boundaries. 



The results of subsections 3.2 and 6.3 tell us that for any large simply-connected 
region R that can be tiled by dominos, height functions associated with random 
tilings of i? will cluster around their average. We furthermore know that this average 
height function depends in a monotone way on the values of the height function 
on the boundary of i?, and is stable under certain kinds of slight perturbations 
of the boundary of R. However, what these theorems do not tell us is whether 
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this dependence is robust under scaling as well. Proposition |l7| tells us that such 
robustness does in fact hold for Aztec diamonds: that is, when one normalizes two 
large Aztec diamonds, one finds that the normalized average height functions are 
very close to one another. That this is true along the boundary is a triviality; 
that it is true in the interior is a much subtler property, known to us only as a 
consequence of Theorem |^. 

We conjecture that scaling-robustness of height functions is true in general. That 
is, suppose i?i,i?2,--. are finite, simply-connected, domino-tileable regions that 
grow without bound, such that suitably rescaled copies of the i?„'s converge to 
some compact subset R* of the plane. Moreover, suppose that the height functions 
associated with the boundaries of the i?„'s, when rescaled by the same respective 
amounts, converge to some function on the boundary of R*. Then we believe that 
the average height functions associated with the R^s, when rescaled, converge on 
the interior of R* as well as on the boundary to some function Ti.. If the boundary 
values behave reasonably (perhaps piecewise smoothness suffices), then Ti should 
be piecewise smooth (with reasonably shaped pieces). 

Under this picture, we view Ti as the solution to a somewhat strange sort of 
Dirichlet problem. We will have more to say about this analogy shortly, but first 
we must leave the issues of large-scale structure (embodied in the average height 
function) and discuss the small-scale structure of random tilings. 

Consider simply-connected regions i?i,i?2,... as above. In each region 
choose a north-going domino space (t„ with normalized location {xmUn) in R*, 
so that {xmUn) {x* ,y*) as n ^ oo, and suppose that the asymptotic renormal- 
ized height function Ti. is differentiable at {x*,y*). Assume that {x*,y*) is in the 
interior of R* and that Ti is "non-extremal" at (x*,y*), in the sense that its tilt 
(s,t) — (^, ^) satisfies \s\ + \t\ < 2. Then we conjecture that the placement 
probabilities at the chosen north-going domino spaces an converge. The arctangent 
formula tells us that the conjecture is in fact true for Aztec diamonds. 

Note that if we were to replace each (T„ by another north-going domino space 
a'^ obtained by shifting it by some fixed vector (i, j) with i + j even, we would get 
the same point {x*,y*) in the normalized limit. Hence, the preceding conjecture 
implies approximate local translation-invariance for the first-order statistics gov- 
erning random tilings of large regions, provided one stays away from the boundary 
(and the tilt is non-extremal). 

This corollary gives us a way to understand the importance of our hypothesis of 
non-extremality. For instance, consider the region shown in Figure it has only 
one tiling, whose local statistics are in no sense governed by any of the statistics 
seen in Aztec diamonds. Taking a suitable limit of such regions one gets a height 
function whose tilt {s,t) satisfies \s\ + \t\ = 2 and hence violates non-extremality. 
Indeed, the statistics do not even exhibit local translation-invariance. (Note also 
that for the Aztec diamond itself, 7i(-, •) is extremal at {x*,y*) if and only if the 
asymptotic entropy at normalized location (x*,y*) is zero, which is the case if and 
only if the asymptotic density of 2 x 2 blocks at normalized location (x*,y*) is 
zero.) 

Having made a conjecture about convergence of first-order statistics, we naturally 
wonder about higher-order statistics as well. We conjecture that in fact all finite- 
order statistics in the vicinity of the points {xn,yn) stabilize as n ^ oo, yielding 
statistics that in some sense "belong" to the limit point {x* , y*) (as long as the tilt at 
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Figure 7. Herringbone pattern. 



{x*,y*) is non-extremal). Then, applying the translation- invariance remark made 
in the preceding paragraph, it follows that each {x* ,y*) determines a process whose 
values are domino tilings of the entire plane. For instance, taking the i?„'s to be 
Aztec diamonds and the point {x* , ) to be the center of the normalized diamond, it 
is natural to conjecture that at the center of the Aztec diamond of order n, the local 
finite-order statistics converge to those of the maximal entropy process mentioned 
in the Introduction. (This special case was conjectured in | JPS |.) Letting {x*,y*) 
vary inside the rescaled temperate zone, we would get a two-parameter family of 
tiling-valued processes; they would all be distinct from one another because they 
would have distinct first-order statistics. The maximal entropy process would be 
unique among these processes not only in having the highest entropy but also in 
being invariant under the full group of lattice translations, rather than merely the 
color-preserving subgroup of index 2. 

It can be shown rigorously that such processes, if they exist, have a combinatorial 
analogue of the "Gibbs property" studied in equilibrium statistical mechanics; that 
is, given a tiling of a cofinite subset of the plane whose finite complement is tileable, 
if one conditions the random process on that particular tiling, then the conditional 
distribution on tilings of the entire plane is uniform. 

Here we leave aside caution and put forward some conjectures about what sort 
of shape the ultimate theory we are striving towards will take. These surmises 
might be false, but we believe they are the natural avenues to pursue in further 
investigations of the theory. 

In the first place, we conjecture that the tiling-valued processes associated with 
the points {x*,y*) will turn out to be ergodic, or indecomposable, in the usual 
sense of the theory of dy namical systems. It is not hard to use the ergodic theorem 
for Z^-actions (see |Kr|) to show that every ergodic, translation-invariant (under 
color-preserving translations), tiling- valued random process determines placement 
probabilities p„, Ps, Pw, and Pe and thence determines a tilt (s,t) — (2(p^ — 
Pe), 2{pn —Ps))- We predict that in those cases where the tilt is non-extremal (i.e., 
|s| + |i| is strictly less than 2), there is in fact a unique ergodic Gibbs measure with 
tilt {s,t). If this were true, it would have many nice consequences; for instance, 
the four numbers Pn,Ps,Pe,Pw would all be determined by the pair (s,f), and thus 
would exhibit only two degrees of freedom, despite the fact that the only obvious 
constraint governing them is Pn+Ps-\-Pe+Pw = 1. A further nice property is that the 
temperate zones of Aztec diamonds would be universal in the sense that they would 
manifest, in the limit, all possible forms of non-extremal local behavior that random 
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tilings of large simply-connected regions can manifest away from boundaries. This 
universality is not peculiar to Aztec diamonds, but instead arises from the fact, 
proved in subsection 6.5, that Aztec diamonds exhibit all possible non-extremal 
tilts. 

An especially nice benefit of the preceding conjecture is that it would open the 
door to a variational approach to the problem of finding the average height function 
H on R* given only its values on the boundary of R*. Given any candidate for H, 
define Nn as the number of domino tilings of i?„ whose normalized height functions 
stay close to Ti. It does not seem too far-fetched to hope that the logarithm of 
Nn, when divided by the area of i?„, converges to an integral over R* , in which 
the integrand is the entropy associated with the unique ergodic Gibbs process with 
tilt (^, Since finding the average height function on i?„ corresponds in some 
sense to maximizing Nn, we would hope that finding the asymptotic normalized 
height function on R* corresponds to maximizing this integral. It might not always 
be possible to solve the associated calculus of variations problem explicitly, but 
such a theorem would be a major advance towards a complete understanding of 
how the presence of boundary conditions can affect the behavior of a domino tiling 
in the interior of a region. 

The preceding idea has in fact been used by physicists, in the context of crystals; 
see for example [NHB, pp. 3562-3563]. There, it is claimed that the shape of a 
crystal surface is determined by minimizing the total surface free energy, which 
is obtained by integrating a local contribution (the surface free energy density) 
depending only on the gradient of the surface. This is believable physically, but 
in any particular lattice model it seems difficult to establish rigorously; it is not 
even clear on purely mathematical grounds why there should exist a surface free 
energy density depending only on the gradient. The analogous statement in random 
tiling theory is the existence of a local entropy depending only on the tilt of the 
height function, but it is conceivable (although we consider it unlikely) that the 
local entropy might not be determined by the local asymptotic behavior of the 
normalized height function. The only approach that we know of that might lead to 
a rigorous proof (or even a heuristic argument) is to prove the conjectures above 
about local statistics and Gibbs measures. 
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